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CHAPTER  1 


INTRODUCTION 

1 . 1  Introduction 

In  recent  years,  interest  has  grown  in  the  development 
and  use  of  automatic  devices  to  test  and  checkout  physical 
systems  of  all  types.  Most  of  the  attention  in  this  field 
has  been  oriented  toward  built-in-test  (BIT)  diagnostic 
subsystems  for  modular  military  electronic  equipment,  mainly 
in  airborne  and  ground  electronic  equipment.  BIT  diagnostics 
have  the  advantage  of  allowing  fewer  and  less  qualified 
maintenance  personnel  and  fewer  pieces  of  external  test  equip¬ 
ment,  which  are  generally  quite  expensive. 

A  primary  equipment  is  composed  of  modular  line 
replaceable  units  (LRUs),  all  of  which  operate  independently. 
Associated  with  each  unit  is  an  a  priori  probability  of  being  in 
failure,  and  it  is  assumed  that  the  probabilities  of  multiple 
failures  are  negligible. 

Whenever  the  equipment  malfunctions,  a  single  LRU  is 
assumed  to  have  failed,  and  two  types  of  diagnostic  tests 
should  be  used  for  the  primary  and  the  secondary  isolations. 

The  primary  isolation  tests  will  be  automatically  executed  by 
the  BIT  in  order  to  identify  the  group  of  LRUs  which  contains 
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the  faulty  unit.  After  the  execution  of  the  automatic  BIT, 
secondary  isolation  will  be  performed  by  semi-automatic  or 
manual  means  which  incur  time  and  extra  equipment  costs  to 
locate  the  single  failed  unit  within  a  group  of  LRUs. 

1 . 2  Statement  of  the  Problem 

Assume  that  equipment  consists  of  n  mutually  exclusive 
groups  of  LRUs.  Associated  with  each  LRU^  is  an  a  priori 
probability  p^,  which  is  the  probability  before  any  diagnosis 
that  the  malfunction  of  the  equipment  is  caused  by  the  failure 
of  LRU^ .  Whenever  the  equipment  fails,  the  BIT  automatically 
executes  a  sequence  of  primary  diagnostic  tests  to  isolate  the 
group  which  contains  the  single  faulty  LRU. 

Each  LRU  could  be  either  good  or  bad,  therefore  a  set 
of  2n  tests  is  required  to  constitute  a  complete  set  of  all 
possible  binary  tests.  However,  if  a  test  that  checks  a 
subset  of  LRUs  is  passed,  the  test  that  checks  the  complement 
of  this  subset  must  be  a  failure,  and  conversely.  Therefore, 
such  a  pair  of  tests  is  redundant,  and  in  the  quest  for  least- 
expected-cost  procedures  the  more  expensive  of  the  pair  can 
be  ignored.  By  this  argument  the  number  of  possible  different 
tests  which  can  exist  for  n  LRUs  is  (2n_^  -  1)  after  excluding 
the  two  tests  which  examine  all  or  none  of  the  n  LRUs. 

Associated  with  each  test,  T^  which  is  included  in  the 
primary  diagnostic,  there  is  a  known  cost,  C^.  The  total  cost 
of  locating  a  particular  faulty  element  is  the  sum  of  the  costs 


of  the  tests  along  the  path  which  leads  from  the  initial  state, 
in  which  no  LRU  is  known  to  be  good  or  bad,  to  the  final 
state  representing  the  group  containing  the  faulty  unit,  plus 
the  cost  of  secondary  isolation  of  that  group. 

The  problem  is  to  design  minimum-expected-cost  test 
procedures  to  be  executed  by  the  BIT.  These  can  be  described 
by  tree  structures,  with  nodes  and  twigs.  Each  node  of  a  tree 
can  be  interpreted  as  a  state  of  ambiguity  subset.  The 
ambiguity  subset  at  each  node  consists  of  the  twigs  that  are 
descendant  from  the  node.  The  test  applied  at  a  node  serves 
to  partition  the  associated  ambiguity  subset,  thus  reducing  the 
ambiguity.  The  root  node,  or  full  subset,  corresponds  to  a 
state  of  complete  ambiguity,  while  at  the  twigs,  which 
correspond  to  unit  subsets  and  hence  where  the  outcome  is 
determined,  there  is  no  further  ambiguity. 

1 . 3  Determination  of  States  Following  a  Test 

The  following  notation  will  be  used  throughout  this 
research. 

T^:  Test  k.  A  test  is  represented  by  an  (n-bit)  number 

containing  only  the  binary  digits  0  and  1.  A  0  is 
assigned  in  position  i  of  a  test  if  LRU^  must  be  good 
in  order  for  the  test  to  pass.  A  1  is  placed  in 
position  i  of  a  test  if  LRU^  is  not  tested. 

/s 

S:  State  of  the  equipment  prior  to  performing  the  test  T^ 

A  state  is  represented  by  an  (n-bit)  number  containing 


n* 


S(S,Tk)  : 


S(S,Tk) : 


n  (S)  : 
C  • 

he* 

T  (S)  : 


T (S,S)  : 


N  (S)  : 
I(S)  : 

b(S,S) : 
(S)  : 


only  the  binary  digits  0  and  1.  The  n  bits  in  the 
designation  of  a  state  ocrrespond,  sequentially  from 
left  to  right,  to  LRU^,  LRl^  •  •  •  LRUn«  A  0  is 
assigned  in  position  i  of  a  state  if  LRU^  is  known  to 
be  good.  A  1  is  assigned  in  position  i  of  a  state  if 
LRUi  is  not  yet  tested.  In  the  initial  state  there  are 
l's  in  all  positions  since  none  of  the  LRUs  have  been 
tested . 

State  of  the  equipment  if  test  Tk  passes.  This 
state  is  computed  by  multiplying  S  and  Tk  bit  by 
bit  with  no  carry. 

State  of  the  equipment  if  test  Tk  fails.  This  state 

A  «. 

is  computed  by  multiplying  S  and  Tk,  the  complement 

of  Tk ,  bit  by  bit  without  carry. 

Number  of  remaining  untested  LRUs  at  state  S. 

Cost  of  test  T,  . 

k 

Set  of  all  possible  tests  which  could  be  used  at 
state  S. 

Set  of  all  tests  which  could  be  used  to  reach  state 
S(S,Tk)  from  state  S. 

A  node  representing  state  S. 

Set  of  indices  of  the  n(S)  remaining  untested  LRUs 
at  node  N(S) . 

A  branch  leading  from  node  N(S)  to  node  N(S). 

Minimum  expected  cost  of  a  sequence  of  tests,  given 
that  the  current  state  is  S. 


C(S,S) :  Expected  cost  of  testing  branch  b(S,S). 


:  '  Expected  cost  for  secondary  isolation  of  LRU^. 


The  basic  structure  of  a  sequential  testing  diagram  is 
illustrated  in  Figure  1.1. 


Figure  1.1.  Sequential  testing  diagram  for  example 
of  4-LRUs . 
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In  order  to  explain  how  the  expected  cost  of  any 
sequential  test  procedure  is  computed,  two  feasible  solutions 
to  the  example  problem  defined  in  Table  1.1  are  shown  in 
Figure  1.2. 


TABLE 

1.1 

EXAMPLE 

PROBLEM 

LRUi 

1 

2 

3 

4 

pi 

.45 

.30 

.20 

.05 

E,- 

6 

3 

5 

1 

Tk 

Binary 

Designation  of 

Test 

Ck 

T1 

1 

0 

0 

0 

$  1 

T2 

0 

1 

0 

0 

6 

T3 

0 

0 

1 

0 

2 

T4 

0 

0 

0 

1 

7 

T5 

1 

1 

0 

0 

3 

T6 

1 

0 

1 

0 

5 

T7  1 

1 

0 

0 

1 

4 

Feasible  Testing  Diagrams  (a)  and  (b) 


Let  be  the  expected  cost  of  the  feasible  test 

procedure  [T1000,  T1100'  T0010^  rePresented  bY  tree  (a)  in 
Figure  1.2. 


Cn  =  C(llll) 


=  C1Q00  +  P1  CdOOO)  +  (p2+p3+p4)  5(0111) 


C1000  +  P1  C(1000)  +  (P2+P3+P4)  (C1100+FvFpT+p7C 


(0100) 


+  (P3+P4) 


P2+P3+P4 


C(0011) ) 


=  C 


1000  +  P1  C(1000)  +  P2  C<0100>  +  (P2+P3+P4} (C1100 


+  (P3~*~P4 ) 


(C 


P2+P3+P4  0010  p3+p4 


C(0010)+ 


4  C(0001))) 


p3+p4 


=  C1000  +  pl  E1  +  P2  E2  +  P3  E3  +  P4  E4 

+  (P2+P3+p4)  C1100  +  (P3+P4)  C0010 

=  7.8 


Let  C2  be  the  expected  cost  of  the  feasible  test 
procedure  [(Tn0(J,  T100Q)  ,  or  (Tnoo,  T001())]  represented  by 
tree  (b)  in  Figure  1.2 


C 


2 


C(llll) 


+  (p.+p,)  C ( 1100)  +  (P,+P.)  C(0011) 


linn 


C1100  +  (pl+p25  (C1000+  p^+p 


-  C(1000)+ 


p2  * 

*  -C(OIOO)) 


p3  * 


Pl+P2 


P4  * 


+  'P3+P4"C0010+  f^100101*  P^'0001” 


3  '*"4 


c  + 
^1100 


pl  E1  +  p2  E2  +  P3  E3  +  P4  E4 


+  (Pl+P2)C1000  +  (P3+P4)C0010 


=  8.9 


This  example  shows  that  the  first  sequential  testing 
procedure  [T10(J0,  ^ll0Q>  T001Q]  is  more  economical  than  the 

second  one  [(T1100'  T1000) '  °r  (T1100'  T0010)J. 


CHAPTER  2 


LITERATURE  REVIEW 

All  works  which  have  been  done  in  the  area  of  the 
optimization  of  fault  detection  and  isolation  procedures  are 
directed  toward  solving  two  basic  problems. 

1.  Generating  a  least  expected  cost  testing  sequence 
to  be  executed  by  the  automatic  BIT  diagnostic  (Primary 
Isolation) . 

2.  Determining  a  troubleshooting  sequence  which 
minimizes  the  expected  cost  of  secondary  isolation  to  locate 
the  single  failed  unit  within  a  group  of  LRUs  identified  by 
the  BIT  primary  diagnostic. 

The  only  other  problem  which  has  been  treated  in  the 
literature  is  the  one  which  restricts  the  repertoire  of  tests 
to  those  which  test  only  single  elements  and  without  even 
assigning  any  probabilities  to  these  tests.  In  this  case,  the 
solution  consists  of  deciding  which  test  to  omit  and  in  what 
sequence  to  perform  the  remaining  tests.  This  problem  can  be 
solved  as  a  machine  setup  problem  as  the  one  in  Glassey  [6] . 

2. 1  Primary  Isolation 

Johnson,  et  al.  [9l  proposed  using  the  information- 
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gain  figure-of -merit  in  order  to  find  a  sequence  of  tests 
that  can  be  executed  by  an  automatic  diagnostic. 

In  spite  of  the  fact  that  this  method  is  easy  to  use, 
it  fails  to  guarantee  optimum  cost  sequence. 

Chang  [  3 ]  used  the  distinguishability  criterion  to 
produce  a  low  expected  cost  testing  sequence  which  is  not 
necessarily  an  optimal  sequence. 

Cohn  and  Ott  [  4 ]  presented  a  recursive  algorithm 
which  is  based  on  the  concept  of  dynamic  programming.  They 
used  set  notation  to  design  a  test  tree.  For  every  possible 
ambiguity  subset,  they  assigned  an  evaluation,  consisting  of 
the  least  expected  cost  of  resolving  that  ambiguity.  The 
evaluation  of  the  subset  of  complete  ignorance  is  the  cost 
of  the  optimal  tree.  This  evaluation  function  is  computed  by 
a  recursion  on  the  number  of  elements  in  the  ambiguity  subsets. 

By  treating  the  equipment  states  as  stages  in  a 
sequential  decision  process,  Sheskin  [11]  applied  probabilistic 
dynamic  programming  to  determine  a  minimum  expected  cost 
testing  diagram.  Using  the  recursive  relationship,  the 
solution  procedure  moves  backwards  stage  by  stage.  The 
solution  procedure  begins  by  equating  the  expected  values  of 
the  terminal  states,  which  corresponds  to  the  groups  into 
which  the  equipment  is  partitioned,  to  the  expected  costs  of 
secondary  isolation  for  these  groups.  At  each  state,  a  set 
of  possible  decisions  consists  of  all  of  the  tests  which  can 
be  performed  is  considered  and  the  optimal  testing  sequence 
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at  this  state  is  found,  until  it  finds  the  optimal  testing 
diagram  when  starting  at  the  initial  state. 

Aly  [1]  constructed  the  problem  as  a  search  tree,  and 
presented  a  branch  and  bound  algorithm  to  find  the  optimal  testing 
sequence . 

2 . 2  Secondary  Isolation 

Gluss  [  7 ]  solved  the  problem  of  having  a  fault 
develop  in  a  system  consisting  of  n  modules  where  each  one 
has  several  elements,  and  that  it  is  required  to  dictate  a 
search  strategy  that  will  optimize  the  search  in  some  fashion 
by  minimizing  a  stipulated  cost  function.  He  developed  a 
model,  which  assumes  that  over-all-tests  of  each  module  may 
be  performed,  and  individual  item  tests  within  modules;  also, 
the  search  is  subject  to  the  constraint  that  before  conducting 
item  tests  the  faulty  module  must  first  be  determined  by 
module  tests. 

Firstman  and  Gluss  [  5 ]  extended  the  work  in  the 
previous  model  of  Gluss,  in  which  the  estimation  of  the 
probabilities  of  faults  lying  in  respective  modules  or 
elements  is  performed  in  a  different  way  from  that  in  Gluss' 
paper:  they  are  computed  from  element  reliability  data  by 

manipulation  of  the  element  failure  rate.  Furthermore, 
consideration  is  given  to  fault  symptoms  that  are  supplied  by 
weighting  the  probabilities  according  to  the  symptoms  infor¬ 
mation. 

All  the  previous  search  models  allow  for  the 
possibility  of  a  test  not  indicating  the  true  state  of  the 
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component  tested.  However  when  Butterworth  [2  ]  tackled  the 
problem,  he  used  tests  that  always  give  the  correct  answer. 

He  developed  several  rules  to  find  the  optimal  sequential 
policies  for  series,  and  parallel  systems  of  independent 
LRUs.  For  a  series  system,  he  indicated  that  the  expected 
cost  for  secondary  isolation  of  the  failed  unit,  given  that 
an  equipment  fault  has  been  isolated  to  this  unit  by  primary 
diagnostic  will  be  minimized  by  removing  and  replacing  the 
LRUs  in  a  nonincreasing  sequence  of  the  values  of  the  ratio 
of  their  probability  of  failure  to  the  average  time  of 
removing  and  replacing  them. 

Butterworth's  rules  fail  to  identify  an  optimal  policy 
for  the  simple  system  where  the  testing  costs  are  identical 
for  all  components.  In  this  case  the  condition  implies  that 
all  the  components  have  the  same  failure  probabilities. 

However,  Halpern  [8  ]  presented  a  simple  adaptive  sequential 
testing  procedure  for  the  k-out-of-n  system  with  equal  cost 
of  all  tests.  This  procedure  covers  the  deficiency  of 
Butterworth's  rules. 

2 . 3  Scope  of  the  Research 

From  the  above  section,  it  is  noticed  that  the  only 
approaches  which  guarantee  an  optimum  testing  sequence  are  the 
recursive  procedure  by  Cohn  and  Ott  [4] ,  the  dynamic  programming 
by  Sheskin  [1],  and  the  branch  and  bound  by  Aly  [1]. 

Capitalizing  on  Aly's  approach,  it  seems  very  promising 
to  formulate  the  problem  as  a  search  tree  and  to  find  the 
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optimal  testing  sequence  using  a  branch  and  bound  approach. 
Using  this  approach  efficiently  could  save  a  lot  of  work  in 
comparison  with  using  the  two  methods  previously  mentioned 
because  of  the  savings  in  the  solution  space  achieved  by 
using  strong  dominance  rules  instead  of  finding  the  optimal 
solutions  among  all  possible  solutions  at  each  node  in  the 
solution  by  dynamic  programming  for  instance,  as  shown  in 
Figure  2.1  for  a  four  LRUs  example  which  uses  many  arcs.  The 
same  problem  could  be  formulated  as  the  search  tree  depicted 
in  Figure  2.2.  However,  all  the  arcs  which  lead  to  any  state 
with  only  one  untested  LRU  and  which  resulted  from  applying 
tests  that  remove  the  ambiguity  of  exactly  one  LRU  are 
omitted  in  order  to  simplify  both  the  network  and  the  tree. 

Also,  by  using  a  good  lower  bound  at  each  node,  most 
of  the  active  nodes  could  be  fathomed  and  the  optimal  solution 
could  be  found  as  fast  as  possible  by  using  a  strong  branching 
rule . 

The  efficiency  of  the  solution  by  using  branch-and- 
bound  approach  depends  upon  the  strength  of  the  bounds,  the 
dominance  and  the  branching  rules.  Consequently,  the  main 
effort  will  be  directed  toward  finding  the  dominance  rules 
which  minimize  the  number  of  branches  as  much  as  possible, 
finding  the  branching  rules  which  concentrate  the  search  only 
in  the  very  promising  branches,  finding  a  lower  bound  at  each 
node  which  helps  in  fathoming  the  maximum  number  of  nodes,  and 
constructing  a  sound  and  efficient  branch-and-bound  algorithm. 
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CHAPTER  3 


A  BRANCH  AND  BOUND  APPROACH 

In  this  chapter  a  branch  and  bound  algorithm  is 
developed  to  find  the  optimal  sequence  of  diagnostic  tests 
to  be  executed  automatically  to  isolate  the  group  of  modules 
which  contains  the  faulty  unit. 

3 . 1  Concept  of  Branch  and  Bound  Techniques 

As  stated  by  Lawler  and  Wood  [10],  branch  and  bound 
is  a  method  of  controlled  search  of  the  space  of  all  feasible 
solutions.  The  space  of  all  feasible  solutions  is  repeatedly 
partitioned  into  smaller  and  smaller  subsets,  and  a  lower 
bound  (in  the  case  of  minimization)  is  calculated  for  the 
value  of  the  objective  function  over  the  solutions  within  the 
subsets.  If  a  known  feasible  solution  is  available  (an  upper 
bound) ,  then  after  each  partitioning  those  subsets  with  a 
lower  bound  exceeding  the  current  upper  bound  are  excluded 
from  further  consideration.  Partitioning  continues  until  a 
feasible  solution  is  found  such  that  its  cost  is  not  greater 
than  the  lower  bound  for  any  subset. 

Branch  and  bound  algorithms  have  two  main  character¬ 
istics;  the  branching  and  bounding  characteristics.  The 
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branching  characteristic  guarantees  that  an  optimal  solution 
will  eventually  be  obtained.  The  bounding  characteristic 
furnishes  the  possibility  of  recognizing  an  optimal  solution 
prior  to  complete  enumeration. 

Therefore,  any  branch  and  bound  algorithm  needs  to 
define  a  set  of  rules  for  (1)  branching  from  nodes  to  new 
nodes,  (2)  determining  lower  bounds  for  the  new  nodes,  (3) 
choosing  an  intermediate  node  from  which  to  branch  next, 

(4)  recognizing  when  a  node  contains  only  infeasible  or  non- 
optimal  solutions,  and  (5)  recognizing  when  a  final  node 
contains  an  optimal  solution. 

In  order  to  use  the  branch  and  bound  technique  to 
find  the  optimal  sequence  of  tests  to  be  used  in  detecting 
and  isolating  the  malfunctioned  unit,  the  search  tree  is 
constructed  as  the  one  in  Figure  2.2  with  a  few  modifications. 
There  is  no  need  for  all  the  nodes  in  the  last  level,  which 
have  states  including  only  one  untested  LRU,  since  their 
status  can  be  found  once  the  search  reaches  any  node  with 
state  of  having  exactly  two  untested  LRUs  regardless  of  its 
level.  Consequently,  savings  can  be  made  in  both  time  and 
storage  required  for  the  solution.  Also,  at  any  node  if  a 
test  which  does  not  remove  the  ambiguity  of  exactly  one  unit 
(in  other  words  it  decreases  the  number  of  untested  LRUs  by  more 
than  one)  is  applied,  the  state  of  this  node  will  be  changed  to 
another  two  states  with  more  than  one  level  difference  between 
them  and  the  given  node.  Therefore,  a  dummy  or  fictitious 


node  will  be  added  after  the  given  node  in  order  to  keep  track 
of  the  two  new  branches  since  the  expected  cost  of  applying  this 
test  should  include  the  expected  costs  of  both  branches. 

The  modified  search  tree  which  is  depicted  in  Figure 
3.1  represents  a  four  LRUs  example  (n  =  4),  with  a  maximum 
number  of  levels  of  (n  -  1).  At  any  state  S  with  n(S)  remaining 
untested  LRUs  there  are  (2n  ~  -1)  branches  emanating  from 

this  state.  Each  branch  represents  a  set  of  possible  tests 
which  could  be  used  to  remove  the  same  ambiguity,  and  conse¬ 
quently  leads  to  the  same  new  states  at  another  level  down  the 
tree.  Since  our  objective  is  to  minimize  the  expected  cost  of 
search  tests,  always,  the  test  with  the  minimum  cost  among  all 
possible  tests  at  every  branch  will  be  considered. 

In  order  to  save  in  the  storage  and  time  requirment, 
which  are  the  main  problems  in  this  kind  of  combinatorial 
problem,  the  nodes  (which  represents  the  states)  of  the  tree 
will  not  be  generated  in  advance,  but  they  will  be  generated 
one  by  one  as  the  algorithm  proceeds.  This  will  not  only  save 
the  number  of  nodes  but  also  the  size  of  information  to  be 
stored  at  each  node.  Once  a  feasible  solution  is  found,  all 
the  remaining  branches  which  are  emanated  from  all  active  nodes 
should  be  checked.  At  the  last  generated  node,  all  the  previous 
branches  and  nodes  which  emanated  from  this  node  and  which  are 
already  examined,  fathomed,  or  had  a  feasible  solution  (which 
is  to  be  stored)  are  to  be  cancelled  and  the  search  is  to 
proceed  in  a  new  active  branch.  By  repeating  this  procedure. 


Level 
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the  number  of  nodes  stored  at  any  time  is  minimum  and  relatively 
very  small. 

In  the  algorithm,  the  search  starts  by  finding  a 
feasible  solution  as  quickly  as  possible  by  moving  directly  down 
the  tree  using  the  branching  rules  from  the  initial  node  at  the 
first  level  to  another  node  in  a  successive  level,  to  a  third 
one  in  a  successive  level,  etc.,  until  finding  a  state  of  having 
only  two  untested  LRUs.  Proceed  upward  in  the  same  branches  in 
order  to  update  the  values  of  the  lower  bounds  at  the  fictitious 
nodes  and,  consequently,  find  the  actual  values  of  the  lower 
bounds  in  the  other  branches  of  all  fictitious  nodes.  This 
procedure  guarantees  finding  a  fast  and  good  feasible  solution 
which  enables  us  to  fathom  efficiently  many  nodes,  especially 
since  the  search  proceeds  in  the  most  promising  branch  at  each 
node  according  to  the  branching  rules  after  applying  the 
dominance  rules,  which  eliminates  as  many  branches  as  possible. 

After  finding  a  feasible  solution,  the  algorithm 
proceeds  by  moving  to  the  last  created  node  and  starting 
branching  and  bounding  as  usual  until  fathoming  all  nodes 
emanating  from  it;  then  going  to  the  second  from  the  last 
created  node  and  so  on  until  fathoming  the  first  node  in  the 
tree.  In  this  case  the  last  solution  cor  ...  sponding  to  the  last 
value  of  the  upper  bound  is  optimal. 
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3. 2  Upper  and  Lower  Bounds 

3.2.1  Upper  Bound  at  the  First  Node 

On  a  minimization  problem — like  the  problem  presented 
here — developing  a  reasonable  initial  upper  bound  on  the 
objective  function  value  is  important  because  it  might  help  in 
fathoming  nodes  before  even  computing  the  first  objective 
function  value  associated  with  a  feasible  search  procedure 
generated  by  the  tree.  In  this  case  the  objective  function 
value  associated  with  any  feasible  procedure  may  serve  as  an 
upper  bound. 

Since  the  maximum  number  of  tests  required  to  find 
the  malfunctioned  unit  among  n  LRUs  is  (n  -  1) ,  using  the 
(n  -  1)  tests  that  have  the  minimum  costs  among  all  tests, 
which  isolate  only  one  LRU  at  the  initial  node,  is  sufficient 
to  find  the  malfunctioned  unit  and  consequently  presents  a 
feasible  search  scheme. 

Noting  that  the  cost  of  the  tests  in  the  objective 
function  should  be  multiplied  by  the  probabilities  of  the 
untested  LRU's  at  each  node  in  order  to  find  the  expected  cost, 
neglecting  the  values  of  these  probabilities  (which  are  less 
than  one) ,  and  taking  into  consideration  the  expected  cost  of 
the  secondary  isolation  of  all  n  LRUs,  results  in  a  value  of  a 
possible  and  reasonable  upper  bound. 

This  initial  upper  bound,  U,  is  defined  as 
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U  =  £  C,  +  E  P.*E.  (3.1) 

ket  k  1=1  1 

where 

Ck  =  cost  of  test  k. 

-  prior  probability  of  failure  of  LRU^. 

=  expected  cost  for  secondary  isolation  of  LRIK. 
t  =  set  of  tests  with  the  minimum  (n  -  1)  costs 

among  the  n  possible  tests  which  isolate  only 
one  LRU  if  they  are  used  at  the  first  node  (i.e., 
tests  with  (n  -  1)  zeros  and  only  a  single  one 
in  the  n  bits  such  as  tests  T1000.  TQ1Q0.  T0010, 
Tqooi  case  of  having  only  four  LRUs) . 

3.2.2  A  Lower  Bound  for  Each  Node 

At  each  node  in  the  tree  a  lower  bound  is  computed  based 
on  the  actual  value  of  the  expected  costs  of  all  tests  used 
prior  to  reaching  this  node,  as  well  as  an  estimate  of  the 
minimum  expected  cost  of  tests  required  to  remove  the  ambiguity 
of  all  the  remaining  untested  LRUs  at  this  node. 

A 

At  any  node  N(S)  applying  test  will  generate  two 

__  A  __  /A  _ 

nodes  N(S)  and  N(S)  corresponding  to  states  S(S,Tk>  and  S(S,Tk) 
respectively  which  arises  two  cases  according  to  the  number  of 
remaining  untested  LRUs  at  each  node. 
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Case  1  Min[n(S),  n(S)]  =  1 

In  this  case,  there  is  no  meaning  of  generating  the 
node  corresponding  to  the  state  of  having  only  one  remaining 
untested  LRU  and,  consequently,  there  is  only  one  branch  to 
be  searched,  assuming  it  is  the  one  starting  with  node  N(S). 
Since  the  remaining  untested  LRUs  at  this  node  are  n(S) 
therefore,  at  most  (n(S)  -1)  tests  could  be  used  to  remove 
their  ambiguity,  and  the  cost  of  these  tests  should  be  multi¬ 
plied  by  the  sum  of  the  probabilities  of  the  untested  LRUs  at 
each  of  the  (n(S)  -  1)  nodes. 

Since  2  _<  n(S)  <_  n,  then  the  minimum  possible  sum  of 
probabilities  to  be  multiplied  by  any  cost  is  the  sum  of  the 
minimum  two  probabilities  of  the  remaining  n(S)  LRUs. 

Let  I ( S )  be  the  set  of  indices  of  the  n(S)  remaining 
untested  LRUs  at  node  N(S),  the  prior  probabilities 
Pn(S)  these  LRUs  are  arranged  in  an  ascending  order  such  that 
<  ~P2 • • • <Pn (g) /  and  C(S)  is  defined  as  a  lower  bound  of 
the  minimum  expected  cost  of  tests  required  to  remove  the 
ambiguity  of  the  n(S)  untested  LRUs  at  node  N(S),  then 

C(S)  =  (p.  +  p,)  •  Z  C. 

-i  r  t  J 

3  cn-i 

where  t  =  set  of  the  n(S)  tests  with  the  minimum  n(S)  costs 
among  all  possible  tests  which  could  be  used  at  this  node. 

A 

Let  T(S,S)  be  the  set  of  all  tests  which  could  be  used 

/V 

to  reach  the  state  of  node  N(S)  from  the  state  of  node  N(S) 
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and  C(S,S)  be  the  minimum  expected  cost  of  the  test  required 

A 

to  reach  the  state  of  node  N(S)  from  that  of  node  N(S),  then 

C(S,S)-  =  min  [CR]  •  I  p,. 

keT (S , S)  iel(S) 

Based  on  the  above  discussion,  a  lower  bound  L(S) 

A 

at  node  N(S)  representing  state  S(S,T^)  could  be  found  by 
computing  the  expression 

MS)  =  L  (S)  -  C  (S)  +  C  (S)  +  C  (S ,  S)  (3.2) 

If  N ( S )  is  the  first  node  in  the  tree  with  state  S 
having  n  untested  LRUs  then, 

n 

L(S)  =  C(S)  +  E  p.  •  E.  (3.3) 

i=l  1  1 

Case  2  Min (n (S) ,n (S) ]  >  1 

A 

In  this  case,  applying  test  at  node  N(S)  will 
generate  two  noaes  which  should  be  both  searched.  Instead,  a 

A  A 

fictitious  node  N(S^)  corresponding  to  a  dummy  state 
will  be  assumed  to  have  resulted  from  applying  test  at 
N (S)  and  will  be  inserted  after  node  N(S).  Then,  the  two 

A 

nodes  will  be  emanated  from  N(S^)  and  generate  the  two 

A  a 

branches  b(S^,S)  and  b(S^,S)  as  shown  in  Figure  3.2. 

A 

Finding  the  lower  bound  at  any  fictitious  node  N(S^) 
is  slightly  different  from  finding  it  at  any  other  node,  since 
it  should  include  C(S)  if  the  search  is  moving  downward  in 


- , 
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branch  b(S£,S),  and  includes  C(S)  if  the  search  is  moving 

A. 

downward  in  branch  b(S£,S).  Arbitrary  in  the  algorithm  the 
search  will  move  first  to  the  node  with  the  maximum  number 
of  remaining  untested  LRUs  among  nodes  N(S)  and  N  (S)  .  Even 
though  both  branches  should  be  searched,  this  procedure  will 
minimize  backtracking  because  it  increases  the  possibility  of 
fathoming  more  branches.  So,  if  n(S)  >  n(S)  the  search  will 

A 

move  downward  in  branch  b(S£,S).  Consequently,  the  lower 

A 

bound  at  the  fictitious  node  N(S£)  could  be  computed  using  the 
expression 

L(Sf)  =  L  (S)  -  C(S)  +  C(S)  +  C(S,Sf)  (3.4) 

Since  there  are  no  tests  required  to  change  the  state 

A  ^ 

of  node  N(S£)  to  the  states  of  nodes  N(S)  and  N(S),  then 
C(S£),  C(S£,S) ,  and  C(S£,S)  are  all  equal  to  zero. 


A  final  word  about  the  lower  bounds.  If  the  search 
reaches  node  N(S)  where  n(S)  =  2,  then  a  feasible  solution 
could  be  obtained.  Let  C  be  the  expected  cost  of  this 
feasible  test  sequence,  then 

C  =  L (S)  -  C(S)  +  min [Cj )  *  I  pf  (3.5) 

jeT(S)  iel(S) 

C  is  the  value  of  the  actual  expected  cost  of  a 
feasible  solution  unless  it  results  from  any  branch  emanated 
from  a  fictitious  node,  in  this  case  it  is  only  a  lower  bound 


of  the  actual  value  of  the  expected  copt  of  a  feasible 
solution.  To  find  the  actual  value  the  search  should  go 
upward  the  tree  to  the  fictitious  node  and  update  its  lower 
bound  by  substituting  the  last  value  of  C  instead  of  C(S)  in 

equation  3.4.  Then,  moving  downward  in  the  other  branch 

/\ 

bfSfjT^),  as  in  Figure  3.2,  until  finding  a  node  with  a  state 
having  only  two  untested  LRUs.  At  this  moment  computing  C 
using  equation  3.5  results  in  the  value  of  the  actual  expected 
cost  of  a  feasible  solution  because  in  this  case  the  cost  of 
the  two  branches  (emanated  from  a  fictitious  node)  has  been  taken 
into  consideration. 

3 . 3  The  Branching  Rule 

The  branching  rule  is  the  criterion  used  at  each  node 
N(S)  to  proceed  the  search  in  one  of  the  possible  (2n^  ^-1) 
branches  where  each  branch  represents  a  set  of  tests  which 
could  be  used  to  change  the  state  of  ambiguity  at  this  node 
to  another  state  in  another  level  down  the  tree.  The  more 
effective  the  branching  rules  are,  the  faster  a  feasible 
solution  could  be  reached  and  consequently  the  less  the  time 
and  speed  required. 

The  branching  rule  used  in  the  branch  and  bound 
algorithm  was  proposed  by  Johnson,  et  al.  [  9 )  as  a  method  for 
constructing  a  good  but  not  necessarily  optimum  sequence  of 
tests  that  can  be  executed  by  an  automatic  diagnostic.  Using 
this  rule  will  improve  the  efficiency  of  the  branch  and  bound 
algorithm  because  it  will  guarantee  finding  a  good  feasible 
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solution  as  fast  as  possible. 

This  rule  uses  the  information-gain  f igure-of-merit , 

F^,  which  is  the  ratio  of  the  ambiguity  removed  by  a  test  T^ 
to  the  test  cost,  C^.  This  rule  is  defined  as  follows: 

a  a 

At  any  node  N(S)  with  a  state  having  n(S)  untested 
LRUs,  by  applying  test  T^,  which  has  a  cost  C^,  either  state 
S(S,Tk)  of  node  N(S)  could  be  reached  if  the  test  passes,  or 

^  A  ^  „ 

state  StSfT^)  of  node  N(S)  will  be  reached  if  it  fails.  Then, 

Fk  -  -  [P  log2  p  +  (1  -  p)  log2  (1  -  p)]/Ck  (3.6) 

where  p  =  I  p  /  I  „  p. 

jel(S)  3  jel(S)  J 

A 

Rank  all  tests  at  node  N(S)  in  a  decreasing  order 
according  to  the  values  of  their  F.  According  to  this  order 
the  tests  will  be  chosen  at  this  node. 

3 . 4  The  Dominance  Rules 

Dominance  rules  could  play  a  very  important  part  in 
determining  the  size  of  the  solution  space  and  consequently 
the  size  of  the  search  tree,  especially  if  it  works  at  the 
root  of  the  tree.  Therefore,  attention  should  be  made  in 
order  to  come  up  with  strong  dominance  rules. 

At  any  node  N(S)  by  applying  test  T^  two  nodes  could 
be  reached;  either  node  of  state  S(S,T^)  or  node  of  state 

_  A  ^ 

S(S,Tk)  with  n(S)  and  n(S)  remaining  untested  LRUs  respectively. 

A  /\ 

Divide  the  set  of  all  possible  tests  T(S)  at  node  N(S)  into 
two  subsets  t,  and  x  where; 


-•  •—  t.t—  ■ 

i  if  A  jam  i'l  d3t  yjyth  i  i  r) 
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T  -  ITklTk  e  T(S),  n  (S)  fl  n  (S)  ] 
and 

T  =  [Tk|Tk  e  T (S) ,  n (S)  =  n(S)] 

Assume  further  that  state  S*  is  the  state  with  the 
minimum  number  of  remaining  untested  LRUs  among  states  S  and 

S. 

{X  pj  ,  if  n  (S)  1  n  (S) 

jel (S* ) 

min[  ^  P j /  S  pj  ]  ,  if  n(S)  =  n  (S) 

jcl(S)  jel (S) 

The  following  theorems  explain  the  dominance  ruler, 
which  will  bo  used  in  the  algorithm.  The  detailed  structural 
proofs  of  all  these  theorems  are  presented  in  Appendix  A. 

Theorem  3.4.1 

At  any  node  N(S),  any  branch  generated  by  a  test  T. 

K 

such  that  Tk  e  t  dominates  any  other  branch  generated  by  a 

test  T  such  that  T  c  t  if: 
m  m 

Ck  =  min[Ci] 
ieT(S) 

and  C.  <  C  *  P(T.  ) 

k  —  m  '  k 

Theorem  3.4.2 

At  any  node  N(S),  any  branch  generated  by  test  Tk  such 

that  T.  e  t  dominates  any  other  branch  generated  by  test  T 
*  in 

such  that  T  e  t  if 
m 
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C.  =  min  iCiJ 
*  ieT(S) 

“0  £k  i  £m  ’  p<Tk> 

Corollary  3.4.1 

Theorem  3.4.2  could  also  be  applied  in  the  opposite 

case,  i.e.,  any  branch  generated  by  test  such  that 

dominates  any  other  branch  generated  by  a  test  Tm  such  that  T^e 

if  Ck  =  min 

ieT(S) 

and  £k  i  £m  '  p|Tk1 
Theorem  3.4.3 

At  any  node  N(S)  with  a  state  having  at  most  four 

remaining  untested  LRUs,  if  test  T^  such  that  Tk  e  t  has  the 

minimum  cost  among  all  tests  which  can  be  used  at  N(S),  then 

the  branch  generated  by  T^  dominates  all  branches  which  are 

generated  by  any  other  test  T  such  that  T  ex. 

*  m  m 

A  summary  of  the  dominance  rules  is  presented  in 

Table  3.1  which  summarizes  the  condition  required  to  make 

a  branch  generated  by  test  T^  at  node  N(S)  dominates  another 

branch  generated  by  test  T  ,  where  C.  =  min  (C . ) . 

m  k  ieT(S)1 


TABLE  3.1 


DOMINANCE  RULES 


! 


T 

m 

Tk 

Tm  e  T 

Tm  e  ? 

Tk  0  T 

£k  5  £m  '  p (Tk> 

'kiV'  <V 

Tk  0  i 

£k  i  £m  •  £’<Tk> 

If  n (S)  <  4 
no  other  condition 
if  required 

If  n (S)  >  4 
no  general  rule 
is  founded 

3 . 5  The  Branch  and  Bound  Algorithm 

In  this  section  the  complete  branch  and  bound  algo¬ 
rithm  for  determining  the  sequence  of  diagnostic  tests  to  be 
executed  automatically  by  the  BIT  to  isolate  the  group  of 
modules  (LRUs)  which  contains  the  faulty  unit  is  given. 

The  input  parameters  are: 
n  =  Total  number  of  LRUs 

T  =  Set  of  all  tests  which  could  be  used 

p^  =  Prior  probability  of  failure  of  LRU^,  i  =  1,2,. 

=  Expected  cost  for  secondary  isolation  of  the 
failed  unit  in  LRU^ 

Ck  =  Cost  associated  with  test  T^ 
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Values  of  the  objective  function  are: 

UB  =  Upper  bound  on  expected  total  cost 
L (S)  =  Lower  bound  on  expected  total  cost  at  state  S 
C  =  Expected  cost  of  a  feasible  test  sequence 
The  parameters  for  creating,  fathoming  nodes  and 
branching  are: 


Step 


ND  =  Current  node  number 

n*  =  Counter  for  nodes  created 

S(S,T^)  =  State  S  generated  by  applying  test  at 

/v 


previous  state  S 

n(S)  =  Number  of  the  remaining  untested  LRUs  at  node  S 

N(S)  =  Node  corresponding  to  state  S 

T(S)  =  Set  of  the  2n^-*-l  possible  tests  at  state  S 

I(S)  =  Set  of  the  n(S)  remaining  untested  LRUs  at 
state  S 

C (S)  =  Lower  bound  of  the  minimum  expected  cost  of 

tests  required  to  remove  the  ambiguity  of  the 
n(S)  untested  LRUs  at  state  S 

Y(S)  =  Set  of  the  remaining  feasible  branches  after 
applying  the  dominance  rules  at  node  N(S) 

(each  branch  could  be  generated  by  at  least 
one  test) . 

£(S)  =  Level  of  node  N(S) 

0  Initialize  the  input  parameter,  let  S  be  the  initial 


state  of  node  N(S),  n(S)  =  n,  MS)  =  1 ,  ND  =  1,  and 
n*  =  1.  Compute  UB  using  equation  3.1  and  L(S)  using 
equation  3.3. 


r 


Step  1 


Apply  a  stopping  test  based  on  the  secondary  isolation 

n  n 


costs.  If  C(S)  +  E  p.  *  E.  >  E  Ei 

i=l  i=l 


Use  the 


secondary  isolation  for  all  LRUs,  stop.  Otherwise, 
go  to  2 . 

Step  2  Use  the  dominance  rules  to  find  Y (S) . 

Step  3  Find  the  information-gain  f igure-of-merit  for  each 

branch  or  test  e  Y(S)  using  equation  3.6.  Rank 
them  in  a  decreasing  order  according  to  the  values  of 
their  F. 

Step  4  Start  branching  using  branch  of  test  with  the 
maximum  F  among  all  tests  in  Y(S)  and  remove  this 
branch  (test)  from  Y(S). 

Step  5  Generate  the  new  two  possible  nodes  by  using  at 

node  S,  denote  them  N(S^)  and  N(S2)«  Find  n(S^)  and 


n(S2) 


Step  6  If  min  [n(S^),  n(S2)]  =  1,  let  node  number  n*  +  1  be 

the  node  with  max  (n(S^),  n(S2)],  go  to  7.  Otherwise, 
let  node  number  n*  +  1  be  a  fictitious  node,  go  to  8. 
Step  7  Let  state  S  be  the  state  of  node  number  n*  +  1,  let 
ND  =  n*  +  1,  MS)  =  MS)  +  1,  go  to  11. 

Step  8  Let  state  S  be  the  state  of  node  number  n*  +  1,  let 

t 

ND  =  n*  +  1,  £ (S)  =  i (S)  +1. 

Find  L(S)  of  the  fictitious  node  N(S)  using  equation  3.4, 
Step  9  If  n(S2)  >  n(S^),  let  node  number  ND  +  1  be  N(S^), 

and  node  number  ND  +  2  be  N(S2).  Otherwise,  let  node 
number  ND  +  1  be  N(S2)  and  node  number  ND  +  2  be  N(S^). 


.1-. 


Step  10 

Step  11 
Step  12 

Step  13 

Step  14 

Step  15 

Step  16 


H (Sx)  =  US)  +1  and  US2)  =  US)  +  1,  let  S  be  the 
state  of  node  number  ND  +  2.  Let  n*  =  ND  +  2. 
Compute  L(S)  using  equation  3.2. 

Apply  the  secondary  isolation  stopping  test.  If 

C (S)  +  Z  pi  '  E-  L  £  Ei  .  Go  to  15.  Other- 
iel(S)  1  iel(S) 

wise,  go  to  13. 

If  L(S)  UB'  fathom  node  N(S).  Go  to  21.  Other¬ 
wise,  generate  T{S),  go  to  14. 

If  n(S)  =  2,  compute  C  using  equation  3.5,  go  to  16. 
Otherwise,  go  to  2. 

Compute  C  =  L(S)  =  C(S)  -  Z  ’  E^  +  Z  E^ , 

iel(S)  iel(S) 

Y (S)  is  empty. 

If  C  >  UB,  fathom  node  N(S),  go  to  21.  Otherwise, 


go  to  17. 


Step  17  If  MS)  =  2,  the  last  solution  is  feasible.  Let 

UB  =  C,  and  state  S  be  the  state  of  node  number  n*, 
go  to  22.  Otherwise,  go  to  18. 

Step  18  If  node  N(S)  is  branched  directly  from  a  fictitious 
node,  go  to  19.  Otherwise,  go  upward  the  same 
branch  to  the  next  node,  let  S  be  the  state  of  this 
node,  with  number  ND. 

Step  19  If  node  number  ND-1  is  fictitious,  let  S  be  its  state 
and  let  ND  =  ND-1,  go  to  17.  Otherwise,  let  S  be 
the  state  of  node  number  ND-2  and  let  ND  =  ND-2,  go 


to  20 . 
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Step  20  Update  the  lower  bound  at  the  fictitious  node  N(S) 
by  substituting  the  last  value  of  C  instead  of 
C(S)  in  equation  3.4.  Let  S  be  state  of  node 
number  ND  +  1,  let  ND  =  ND+  1.  Compute  L(S),  go  to 
12. 

Step  21  Let  ND  be  the  number  of  the  node  N(S).  If  N(S)  is 
branched  from  a  fictitious  node,  fathom  also  node 
number  ND-1 ,  and  let  S  be  the  state  of  node  number 
(ND-2 )  and  let  its  number  be  ND.  Otherwise,  let  S 
be  the  state  of  node  number  (ND-1)  and  let  its 
number  be  ND. 

Step  22  If  MS)  =  1,  go  to  30.  Otherwise  go  to  23. 

Step  23  If  MS)  =  2,  go  to  28.  Otherwise  go  to  24. 

Step  24  If  N(S)  is  fictitious,  or  n(S)  =  2,  let  state  S  be 
the  state  of  node  number  ND-1,  and  its  number  is 
ND,  go  to  22.  Otherwise  go  to  25. 

Step  25  If  Y(S)  is  empty,  go  to  26.  Otherwise  go  to  31. 

Step  26  If  N (S)  is  branched  from  a  fictitious  node,  go  to  27. 

Otherwise,  let  state  S  be  the  state  of  node  number 


ND-1  and  let  its  number  be  ND,  go  to  22. 

Step  27  If  the  lower  bound  at  the  fictitious  node  has  been 
previously  updated,  let  state  S  be  state  of  node 
number  ND-1  and  let  its  number  be  ND,  go  to  22. 
Otherwise,  go  upward  this  branch  to  the  next  node 
let  its  state  be  S  and  its  number  ND,  go  to  22. 
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Step  28  If  N(S)  is  fictitious,  or  n(S)  =  2.  Let  S  be  the 
initial  state,  go  to  3-.  Otherwise  go  to  29. 

Step  29  If  Y ( S)  is  empty.  Let  S  be  the  initial  state,  go  to 

30.  Otherwise  go  to  31. 

Step  30  If  Y ( S )  is  empty,  stop,  go  to  32.  Otherwise  go  to 

31. 

Step  31  Let  n*  be  the  number  of  node  N(S),  go  to  4. 

Step  32  The  optimal  sequence  of  tests  is  the  one  associated 
with  the  last  value  of  the  upper  bound  UB. 

3 . 6  Verification  of  the  Algorithm 

The  algorithm  of  Section  3.5  was  coded  in  FORTRAN  IV. 
The  code  was  verified  using  the  example  problem  used  in  [11] 
and  shown  in  Table  1.1. 

The  search  tree  used  in  solving  this  problem  by  branch 
and  bound  algorithm  is  presented  in  Figure  3.3.  The  same 
optimal  sequence  of  tests  has  been  obtained.  Either  sequence 

of  tests  T1000,  T110Q,  and  TQ010  or  ^10Q0>  TooiO'  and  T1100 
produced  the  same  optimal  solution. 

From  the  tree  presented  in  Figure  3.3  it  is  noticed 
that  the  dominance  rules  and  lower  bounds  worked  efficiently 
to  reduce  the  size  of  the  tree  to  include  only  seven  nodes 
compared  with  the  original  possible  tree  for  four  LRUs,  which 
has  23  nodes  as  well  as  the  dynamic  programming  network  which 


has  16  nodes. 
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It  is  noticed,  also  that  even  though  the  total  number 
of  nodes  generated  during  the  execution  of  the  algorithm  was 
seven,  the  maximum  number  of  nodes  stored  at  any  time  was 
only  three,  which  is  relatively  small  and  reasonable.  Also, 
the  first  feasible  solution  happened  to  be  the  optimal  solution 
which  shows  the  strength  of  the  branching  rules  and  its 
effectiveness  in  helping  fathoming  the  remaining  active  nodes. 

Thus,  in  a  simple  test  example,  the  efficiency  of  the 
algorithm  was  verified. 

3.7  The  Heuristic  Algorithm 

The  branch  and  bound  algorithm  explained  in  section 
3.5  finds  efficiently  the  optimal  solution.  However,  the  size 
of  the  problems  which  could  be  solved  by  this  algorithm  is 
relatively  small  because  of  the  storage  burden  and  time 
requirement,  which  is  inherented  in  most  combinatorial 
problems. 

This  heuristic  algorithm  is  simply  the  same  branch 
and  bound  algorithm  explained  in  the  previous  section  with 
two  more  stopping  tests  which  stop  the  search  for  optimality 
by  stopping  the  search  either  directly  after  finding  the 
second  feasible  solution  or  after  generating  a  limited  number 
of  nodes  based  on  the  maximum  number  of  nodes  required  to 
find  a  feasible  solution.  By  experiment  it  was  found  that 
the  best  results  happened  when  the  search  stopped  after 
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generating  a  number  of  nodes  equals  to  fifteen  tiroes  the 
maximum  number  of  nodes  required  to  find  a  feasible  solution. 
These  tests  were  developed  from  the  computational  results  of 
the  branch  and  bound  algorithm  which  showed  that  most  of  the 
search  time  was  consumed  in  proving  optimality  not  in  finding 
the  optimal  solution  itself. 

The  stopping  test  based  on  the  second  feasible  solution 
can  be  added  in  step  17  in  the  branch  and  bound  algorithm. 

While  the  stopping  test  based  on  the  total  number  of  nodes 
could  be  added  before  step  11. 

The  value  of  the  objective  function  obtained  by  the 
heuristic  algorithm  was  found  to  be  on  the  average,  99.244% 
or  more  of  the  values  of  the  optimal  solution  for  all  test 
problems.  The  details  of  the  computational  experience  are 
presented  in  Chapter  4. 


CHAPTER  4 


COMPUTATIONAL  RESULTS 

In  this  chapter  the  computational  experience  with 
both  the  branch  and  bound  and  heuristic  algorithms  presented 
in  Chapter  3  is  demonstrated  and  analyzed.  The  test  problems 
were  randomly  generated  from  uniform  distribution.  All 
probabilities  of  failure  of  the  n  LRUs  were  generated  from  a 
uniform  (0-1)  distribution.  The  costs  of  all  tests  were  generated 
from  a  uniform  (1-20) ,  while  the  expected  costs  for  secondary 
isolation  of  all  LRUs  were  generated  from  a  uniform  (1-10) 
distribution.  All  problems  were  run  on  the  University  of  Oklahoma 
IBM  370/158J  computer.  The  results  are  summarized  in  Table  4.1. 

As  in  all  combinatorial  problems,  the  required  compu¬ 
tational  time  is  a  function  of  the  size  of  the  problem  as  well 
as  the  number  of  active  nodes.  As  depicted  in  Figure  4.1  the 
case  of  n  >  8  LRUs  is  the  critical  case  where  the  time  starts 
increasing  exponentionally  from  10,6  seconds  in  case  of  n  =  7 
to  160.345  seconds  in  case  of  n  =  8. 

Table  4.2  displays  a  comparison  between  the  branch  and 
bound  algorithm  and  the  heuristic  one.  The  savings  in  compu¬ 
tation  time  by  using  the  heuristic  algorithm  is  obvious, 
especially  when  the  number  of  LRUs  increases.  However,  the 
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TABLE  4.2 

A  COMPARISON  BETWEEN  THE  RESULTS  OF  BRANCH  AND 
BOUND  AND  THE  HEURISTIC  ALGORITHMS 


LRU 

Average  CPU  Time 
Sec 

%ge  of  time 
saved  by 
using 
heuristic 

%ge  of  difference  in  the 
objective  function  between 
the  optimal  and  heuristic 

B&B 

Heuristic 

Average 

Maximum 

4 

.331 

.318 

3.86% 

.4418% 

3.91% 

5 

.526 

.426 

18.95% 

.756% 

7.25% 

6 

1.765 

1.032 

42.17% 

.349% 

3.37% 

7 

10.6 

6.113 

42.45% 

.  666% 

3.4% 

8 

160.345 

23.579 

85.3% 

.0124% 

.0749% 

9 

1253.715 

74.137 

94.00% 

0 

0 

10 

>3600 

280.135 

>92 . 2% 

optimal  s 

solution 

is  not  known 

sacrifice  in  the  optimal  value  of  the  objective  function  is 
less  than  7.25%  of  the  optimal,  and  on  the  average  it  is  less 
than  0.756%.  Also,  Figure  4.2  shows  that  the  heuristic 
algorithm  reached  the  optimal  solution  in  more  than  82%  of  the 
problems  tested  which  shows  the  effectiveness  of  this 
algorithm. 

Table  4.3  shows  a  comparison  between  the  branch  and 
bound  algorithm  and  the  dynamic  programming  approach  used  in  [11]. 
This  comparison  is  based  on  the  maximum  number  of  nodes  created  by 
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TABLE  4.3 

A  COMPARISON  BETWEEN  THE  BRANCH  AND  BOUND 
AND  DYNAMIC  PROGRAMMING  ALGORITHMS 


V  LRU 

Mr 

m.  oi 
nodes 

created  N, 

3 

4 

5 

6 

7 

8 

9 

Branch  and 

Bound 

4 

26 

113 

432 

666 

2812 

1959 

Dynamic 

Programming 

i _ 

12 

77 

39 

1767 

7560 

31369 

128010 

using  both  algorithms.  The  comparison  indicated  a  dramatic 
difference  in  the  number  of  nodes  created  especially  for 
n  ^  7  LRUs.  A  comparison  in  the  computation  time  would  have 
been  rather  more  important.  However,  no  computational  results 
were  reported  in  case  of  using  dynamic  programming,  only  the 
upper  bound  on  the  number  of  states  generated  by  dynamic 
programming . 


CHAPTER  5 


SUMMARY  AND  CONCLUSIONS 


5 . 1  Summary 

Two  approaches  to  the  cost  effective  design  of  fault 
isolation  procedures  were  investigated.  The  problem  was 
formulated  as  a  search  tree  in  which  the  optimal  search 
procedure  could  be  found  using  a  branch  and  bound  approach. 
Dominance  and  branching  rules  were  developed,  then  a  branch  and 
bound  algorithm  was  presented. 

Having  studied  the  computational  results,  another 
heuristic  algorithm  was  developed  which  proved  to  be  efficient 
and  fast.  An  example  problem  was  solved  to  illustrate  the 
efficiency  of  the  branch  and  bound  algorithm  and  was  compared 
with  a  previous  dynamic  programming  algorithm. 

Computational  results  indicated  that  the  heuristic 
algorithm  was  faster  than  the  branch  and  bound  one  with  a  very 
slim  sacrifice  in  optimality. 

Computational  results  were  reported  and  compared  to  the 
available  results  of  other  algorithms. 

5.2  Conclusions 


Several  conclusions  can  be  drawn  from  this  research 
regarding  the  consideration  of  new  approaches  for  fault  isolation 
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problems.  They  are: 

1.  The  branch  and  bound  approach  could  be  used 
successfully  to  tackle  the  problem  of  designing  a  cost 
effective  fault  isolation  procedure.  Because  of  the  branching 
and  dominance  rules,  many  of  the  nonoptimal  solutions  would  be 
eliminated  early  in  the  solution  procedure  which  could 
efficiently  reduce  the  size  of  the  required  search  tree,  as 
well  as  the  time  and  storage  needed  to  find  the  optimal 
solution. 

The  branch  and  bound  algorithm  proved  to  be  more 
efficient  than  the  dynamic  programming  scheme  which  has  been 
used  in  previous  works  to  seek  optimal  procedures. 

3.  The  heuristic  algorithm  presented  in  section  3.7 
proved  to  be  a  good  compromise  between  the  ultimate  goal  of 
optimality  and  the  problem  of  time  requirement  to  achieve  this 
goal.  This  algorithm  has  the  advantage  of  finding  a  near 
optimal  solution  in  a  very  short  time  compared  to  other 
methods . 

4.  Even  though  the  size  of  problems  solved  efficiently 
by  the  two  algorithms  are  limited  to  nine  LRUs,  this  size  is 
still  greater  than  any  problem  reported  to  be  solved  in  any 
previous  work. 

5 . 3  Future  Work 

Recommendations  for  further  research  in  the  cost 
effective  design  of  fault  isolation  procedures  would  be: 


_  — ~ 
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1.  Developing  a  technique  to  minimize  and  control 
the  number  of  possible  tests  which  could  be  used  in  the 
search  because  of  the  dramatic  increase  of  the  possible 
number  of  tests  with  the  increase  of  LRUs. 

2.  More  investigation  in  developing  branching  and 
dominance  rules  and  more  work  in  designing  test  procedures 
using  branch  and  bound  approaches. 

3.  Investigating  how  to  partition  the  equipment  into 
optimum  groups  of  modules. 

4.  Considering  the  problem  without  neglecting  the 
possibility  of  multiple  failures  of  two  or  more  LRUs  at  the 
same  time. 

5.  Studying  the  effect  of  imperfect  information  on 
the  optimum  test  procedures  and  how  to  modify  the  solution 
according  to  that  (sensitivity  analysis) . 

6.  Determining  an  optimum  procedure  which  minimizes 
the  expected  cost  of  secondary  isolation  to  locate  the  single 
failed  unit  within  the  group  of  LRUs  identified  by  the  BIT 
primary  diagnostic. 
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APPENDIX  A 


DOMINANCE  RULES  THEOREMS 


The  proofs  of  all  theorems  which  have  been  used  to 
determine  the  dominance  rules  in  Section  3.4  are  presented 
here . 

In  order  to  simplify  the  proofs,  the  cost  of  a 

test  T^  which  has  l's  in  positions  i,j,...,z  will  also  be 

identified  as  C .  . 

i /3/***/Z 


Theorem  3,4.1 

At  any  node  N(S),  any  branch  generated  by  a  test  T^ 

such  that  Tj^  c  t  dominates  any  other  branch  generated  by  a 

test  T  such  that  T  e  t  if: 
m  m 

Ck  =  minted] 
ieT(S) 

and  '  P(V 

Proof 

Referring  to  Figure  A.l  and  Figure  A. 2  which  represent 
two  branches  from  the  search  tree  of  a  problem  of  5  LRUs,  both 


52 


Figure  A.l  A  Search  Tree  for  a  5-LRUs  Examp 


branches  are  emanated  from  the  first  node  and  by  using  tests 
belong  to  set  t.  All  tests  which  could  be  used  at  any  branch 
are  presented  in  Table  A.l. 

Let  the  cost  of  branch  b^(SQ,S^,S^)  which  pass  through 
nodes  N(SQ),  NfS-^)  and  N(Sg)  be  Cb  (S^S^Sg)  ,  and  let  the 
cost  of  branch  b2  (S^ , Sg , , S^g)  which  pass  through  nodes  N(Sq), 
N(S6),  N(Sy),  N(S1q)  be  Cfc  (SQ , Sg , S? , S1Q) ,  then, 

Cb  (SQ,S^,Sg)  C4  g^"(p4tpg)  *  min  [C^  ,C4  ,C^  4f^3  4*^1  4* 

C2,5'C3,5J+(pl+P2+P3)  *  min[Cl'Cl,4'C2,3'Cl,5J+ 
(p2+p3)  -gmintC^  2,C2,C3,C2  4»C2  g,C3  2,C3  4' 

C  ]  +  E  p .  E  . 

i=l  1  1 


Cb2(S0'S6'S7'S10)  =  C1+(P2+P3+P4+P5)  *  min [C1 , 2 'C2 1  +  <P3+P4+P5} 

•  min[C2  3/clf 3'c3'c4f 5J+(P4+P5)  *  rointCg, 

C4 ,C1 , 4 ,C3 , 4 ,C1 , 5'C2 , 4 ,C2 , 5'C3 , 5^ +if 1Pi  *  Ei‘ 


Branch  b^(Sg,S3,Sg)  generated  by  test  T 


00011 


dominates  branch 


b2  ^S0'S6'E7'S10^  9enerated  bY  test  T10000  if 

C4  5  +  (pl+p2+p3^  ’  +  (P2+P3 )  '  minlCj^  2,C2, . . .  ] 

<  C1  +  (p2+p3+p4+p5)  •  min[C1^2,C2l  +  (p3+p4+p5)  *  min[C4  5,...] 

But  if  C.  ;  =  min  [C . ] 
ieT (SQ)  1 


and  since  E  p.  =1,  therefore  ... 
i=l  1 


•f 
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( 


branch  b^iS^S^Sg)  dominates  branch  b2  (SQ,Sg,S7 ,S1Q)  if 

-  (P4+P5)C1  <  (P4+P5)  ’  min[C1>2,C2]  -  (P1+P2)c4f5 

i.e>  (P4tpg)C^  ^  (P^+P2 ^4 f  5  —  (Pjj^Pg)  ’  min  [C^  ^  2  < ^2 ^ 

which  could  be  satisfied  if  C,  c  <  (P«+pc)C, . 

4,5  —  r4  5  1 


This  proof  is  valid  even  if  the  cost  of  branches 
^(SqjS^S^)  and  b4  (Sq,S^,S^)  are  less  than  that  of  b^(Sg, 
Sg,Sg)  because  in  this  case  they  dominate  branch  b^  (Sq,S^  ,Sc_) 
and  consequently  dominate  branch  b2  (Sq , Sg , S7 , S^q) .  However, 
it  is  not  the  same  for  other  branches  bg (Sq , Sg , S7 , S32)  and 
b2 (Sg,Sg,S7 ,S^q) .  Therefore,  it  should  be  proved  that  branch 
b^  (S^S^Sg)  dominates  branch  b2 (SQ , Sg , S7 , S1Q)  and  any  other 
branch  generated  by  test  T^ooOO'  whatever  the  branch  emanating 
from  node  N(Sg)  with  minimum  cost  is. 


Case  1  If  branch  bg (Sq , Sg , S7 ,Sg2)  is  the  optimal  branch 
generated  by  test  T^0000 

Let  the  cost  of  branch  bg (SQ,Sg,S7 ,S^2)  be  (S^,Sg, 


S7'S12) 

Cbg(S0'S6'S7'S12) 


C1  +  (p2+p3+p4+p5)  •  min[C1>2,C2]  +  (p3+p4+ 

Pg) •min[Cg,C2^g,C1>g,C3^4]  +  (?3+P4) -min [Cg , 


'C,  A 


-C, 


'41, 31, 42, 32, 4'3, 54, 5  ,I,-i  ~i 


.]  +  Ip. 

i=l 


So,  if  C.  ,  <  min[C.] 
isT (Sq) 

then  branch  b^(Sg,Sg,Sg 
C4 , 5  +  (P4+P5)  *  min[Cg 


)  dominates  branch  bg  (Sq , Sg , S7 , Sg2 
'C4'C1,4'C3,4'C1,5'C2,4'C2,5,C3,51 


)  if 
+ 
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(pl+P2+P3)  ’  n'intCl'Cl,4'C2,3'Cl,51  +  (p2+P3)  ’  min  Ici  f  2  ,C2  ,C3  ' 
C2,4'C2,5'C1,3'C3,4'C3,5J  -  C1  +  (P2+P3+P4+P5}  ’  mln tCl f 2 'C2J 
+  (P3+P4+P5)  *  min^C5'C2,5'Cl,5,C3,41  +  (P3+P4}  C4 , 5 
or  -  (p4+p5)C1  <  -  (P1+P2+P5)c4,5  +  (P4+P5)  '  rnin  lc1  /  2,C2] 

+  p3  •  min [C5 ,  . . . ] 

or  (p4+p5)C1  >  (P1+P2+P5)c4,5  "  (P4+P5)  *  It'inICi,2'  •••! 

-  p^  •  min [C5 ,  . . . ] 

which  could  be  satisfied  if 

C4,5  i  IP4+p5>Cl 

under  the  condition  that 

C4  5  =  min[Ci] 
ieT {SQ) 

By  the  same  procedure  it  could  be  proved  that  branch 
b^(Sg,S^,Sg)  dominates  branch  bg  (SQ,Sg,S7  ,S^)  and  any  other 
branch  generated  by  test  Ti0OOO  an<3  a  test  such  that  kex  at 
node  Sg. 

Case  2  If  branch  b^ (Sq , Sg , Sg , S^)  is  the  optimal  branch 
generated  by  T10000 

Lfit  ths  cost  of  brsnch  (SQ^S^^S^/Sg)  (Sq^S^/Sq^s^ 

Cb7 (S0,S6,S8'S9)  =  C1  +  (p2+p3+p4+p5)  '  min  lC2,4'C3,51  + 

(p2+p4 )  ’  rnin  [C^  +  (p3+p5)  •  min 

tC4, 5' * ' - 1 

So  if  C .  =  min  [C .  ] 

ieT (SQ) 
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branch  b1(SQ,S1,S5)  dominates  branch  b7  (Sg,Sg,Sg,Sg) 
if  -  (P4+P5)C1  <  -  px  C4>5 
i.e.  (p4+p5)C1  >  p1  C4^5 
which  could  be  satisfied  if 
C4 , 5  i  <P4+P5)C1 

Therefore,  in  any  event  the  branch  generated  by  Tqqq-^ 
such  that  Tgoon  e  T  dominates  any  other  branch  which  is 
generated  by  T1000()  such  that  T1000Q  e  t  if 

CA  c  =  min  [Cj] 
ieT(S07 

and  C4^5  <  Cx  •  P(T00011) 

Theorem  3.4.2 

At  any  node  N(S),  any  branch  generated  by  test  T^  such 

that  T^  £  x  dominates  any  other  branch  generated  by  test  Tm 

such  that  T  e  t  if 
m 

c,  =  min  iCil 
*  ieT(S) 

and  Ck  <  Cm  -  P(Tk) 

Proof 

Referring  to  Figures  A. 3  and  A. 4  which  represent  two 
branches  in  a  search  tree  of  problems  having  6-LRUs,  the  first 

branch  in  Figure  A. 3  is  generated  by  T^qoOOO  where  T100000  e  T 

and  the  second  branch  in  Figure  A.  4  is  generated  by 

where  £  x.  All  tests  which  could  be  used  at  any  branch 


1 


S) N  f  TITII 
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are  presented  in  Table  A. 2. 

Let  cost  of  branch  b^ (Sq , , S3 , S^)  be  (Sq,S3,S3,S^) 

Cb  (S(),S1,S3,S5)  =  +  (p2+p3+p4+p5+p6)  *  minfC2,3'Ci,2,3^ 

+  (p2+P3^  '  min [C2 »C3 , . . . ,C^  3  g]+cost  of  optimal 

6 

search  starting  from  node  N{S,)+  I  p.  •  E.  . 

J  i=l  1 

Let  cost  of  branch  b2 (SQ,Sg,Sg)  be  (SQ,Sg,Sg) 

Cb2(S0'S6'S8}  =  C1 , 2 , 3  +  (P1+P2+P3)  *  min[Cl'  +  (P2+P3)  * 

min [C_ ,C, , . . . ,C,  .  ,]  +  cost  of  optimal  search 
2  3  L  ,  3  ,  b  g 

starting  from  node  N(S7)+  E  p.  •  E. . 

'  i=l  1 

So  if  =  min[C^] 
ieT(S0) 

then  branch  b1  (S^S^S^Sg)  which  is  generated  by  test  T100000 
where  Ti00000  e  T  dominates  branch  b2  (SQ,S^,SQ)  which  is 
generated  by  test  T^^^000  wbere  TmooO  e  T  ^ 

C1  +  (P2+P3+P4+P5+P6)  *  min  (C2,3'  Cl,2f3)  -Cl,2,3  + 

(pl+p2+p3 )  *  min  ICj,  ...1 

or  if  CL  +  C1>2>3  (p2+p3+p4+p5+p6)  <  C1/2/3  +  <VP2+P3)C1 

-P1  C1 , 2 , 3  -  "  (p4+P5+P6)Cl 

px  c1^2  3  >  (p4+p5+p6)  •  c1 

which  is  still  satisfied  if 


C1  -  P1  C1 , 2 , 3 
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This  proof  is  valid  whatever  the  optimal  branch  from 


node  N(S^)  is  because  if  any  other  branch  dominates  branch 
b1(SQ,S1,S3,S5) ,  it  will  also  dominate  branch  b2 (SQ,Sg,Sg) . 
However,  to  check  the  validity  of  the  proof  if  any  other  branch 
from  node  N(Sg)  is  optimal,  let  the  cost  of  branch  b3(Sg,Sg,Sg) 
be  Cfa3  (SQ,S6,S9) . 

Cb  (Sq/S^/S^)  =  Cx  2  3  +  (P^+P2+P3 )  *  [C^/C^  2/^3  4*^3  5/^3  6 

Cl,2,6'Cl,2,4'Cl,2]+(pl+P2)*min[Cl'**']+  COSt  of 

6 

optimal  search  starting  from  node  N(S7)+  I  p.*E. 

i=l  1 

So  if  C,  =  mintC^] 
ieT (SQ) 


Then  branch  b^ (SQ , , S^)  dominates  branch  b3(SQ, 
Sg,S9)  if  -  p^  2  3  +  (P2+P3^  ‘  min  2'  ^3  4 •  ^3  5* 

c3  g,^  2  g,^  2  4rC1  25/  ...]  ±  (p1+p2)C1  +  (p^+P2+P3)  *  min 

[C3 ,Ci , 2 /C3 , ^ ,C3 , 5 ,C3 , g ,Ci , 2 , & ' ^1 ^  2 , 4 '^1 , 2 , 5^ 

or  p^  C3  2  3  —  *  ^ P3 "^"P4^"5"^"P 6 ^  ^1  P1  *  min  [C3 /  •  •  •  3 

or  P1  Cl,2,3  -  (P3+P4+P5+P6)C1  "  P1  *  min  [C3' 
which  could  be  satisfied  if 

C1  -  P1  Cl,2,3 

By  using  the  same  procedure,  it  could  be  concluded 
that  a  branch  generated  by  test  T^qooOO  where  Ti00000  £  T 
dominates  any  other  branch  which  is  generated  by  test  TmQ00 
where  Tinooo  c  t  if 
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and 


C1  =  min  [C.J 
ieT (SQ) 

C1  -  Cl,2,3  ’  p(T100000) * 


Corollary  3.4.1 

Theorem  3.4.2  could  also  be  applied  in  the  opposite 

case,  i.e.,  any  branch  generated  by  test  T^  such  that  T^ex 

dominates  any  other  branch  generated  by  a  test  Tm  such  that  Tmex 

if  C,  =  min  [C .  ] 
k  l 

ieT(S) 

and  Sk  -  Sn  '  p[V 


Proof 

Referring  to  the  proof  of  theorem  3.4.2,  and  Figures 


3.5  and  3.6 

if  Cl,2,3  i  min  l5i> 
ieT (SQ) 

Branch  b2(S0,Sg,Sg)  which  is  generated  by  test 
where  TmQ00  e  *  dominates  branch  b^  (Sg ,  ,  Sg  , Sg)  which  is 

generated  by  test  T100000  where  T100000  e  T  if 

-(p4+p5+p6)  Cj  <  -  Plcli2i3 

or  (p4+p5+p6)C1  >  PlC1>2i3 


which  could  be  satisfied  if  2  3  i  (P^j+Pg+Pg)*^  • 

This  proof  is  valid  whatever  the  optimal  branch  from 
node  N(Sg)  is.  However,  to  check  the  validity  of  the  proof 
if  any  other  branch  from  node  N(S^)  is  optimal,  let  the  cost  of 
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branch  b4  (Sq/S^S^S^  be  (so'Sl'S2'S4^ ’ 

6  6 

Cb4(S0'Sl'S2'S4>  =  C1  +  Ll2  Pi  ’  min[C2'Cl,2]  +  .f3  pi  *  minlC3' 

Ci  3/^^  2  3'C2  3 1  t crost  of  the  optimal  search 

starting  from  node  having  the  state  000111]  + 

6 

.  r  p .  •  e 

i=lFi  *i- 

So  if  Ci  2  3  =  min  lCj_] 
ieT (SQ) 


Branch  b2(SQ,S6,Sg)  dominates  branch  b4 (SQ, Si , S2 , S4) 
if  Ci  2  3  t  E  *  min  [  C 1 ,  *  •  •  ]  (P2*^"P3^  min  [C2,Ci^2i 


i=l 


6  6 

<  C,  +  E  p.  •  min  [CofC,  0]  +  E  p-  ’  min  [C,  7  , ,  ...] 

—  1  .  j  1  2  1/2  i_2  1  -Lf2,J 

or  -  (P4+P5+P6)  •  Cj  <  -  (P1+P2+P3)C1,2,3 

(P4+P5+Pg)C1  >  <Pl+P2'fP3^Cl,2,3 

which  could  be  satisfied  if 
Cl,2,3  i  <P4+P5+P6>C1 
which  is  still  satisfied  if 

Cl,2,3  -min  UPi+P2+P3),  (P4+P5+P6)l  ‘  Ci 

Therefore,  in  any  event  a  branch  generated  by  test 

T111000  where  TlllOOO  £  *  dominates  any  branch  which  is 
generated  by  test  T100000  where  Ti00()00  e  x  if 

C1 , 2 , 3  =  min  [5i] 
ieT (SQ) 

and 

C1 , 2 , 3  -  C1  *  p (T111000) 
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Theorem  3.4.3 

At  any  node  N(S)  with  a  state  having  at  most  four 

remaining  untested  LRUs,  if  test  T^  such  that  T^  e  x  has  the 

minimum  cost  among  all  tests  which  can  be  used  at  N(S),  then 

the  branch  generated  by  T^  dominates  all  branches  which  are 

generated  by  any  other  test  T„  such  that  T  ex. 

m  m 


Proof 

With  reference  to  Figure  A. 5,  depicting  a  search  tree 
for  a  four  LRUs  example,  all  tests  which  could  be  used  at  any 
node  in  this  tree  are  presented  in  Table  A. 3 

Br anche s  b^(SQ,S^,S2)  ,  ^2^E0,E3f^4^  1  ^3 

emanated  from  the  first  node  N(Sg)  by  tests  T^q^,  T1001 

respectively,  let  the  costs  of  these  branches  be  (Sg,S^,S2), 

C.  (Sn , S,  ,  S . )  ,  and  C,  (SA,S C,S,)  respectively  where  all  the 

U  J  4  D  ^  U  3  o 

tests  belong  to  set  x 
4 

Cb  (Sq , , S2> =  1  P^  *  +  Cx  ^  ^  (P2+P4 )  min [C^ »C2 ,C^  ^ , 

cl,2)  +  (Px+P3)  •  minlC1,C3,Clf4,C12] 

4 

Cfe  (S0,S3,S4)  =  I  Pi  '  Ei  +  ci  2  +  (Pi+P2 )  '  min[C1,C2, 

C1  3'C1  4^  +  (p3+p4)  *  min  [C4  ,C3,C1  3(C^4] 


C,  (S 
b3 


0'S5,S6)  “  if1  pi  *  Ei  +  C1 , 4  +  (P2+P3J  "  mintc3'C2' 
C1,2'C1,3]  +  (P1+P4)  ‘  minfc1'C4'c1>2'Cl,3J 


So,  if  C,  ,  =  min  [C^] ,  then 
i£T(SQ) 


1111  )  N (S 


1010 

T1000'  T0010'  T1001'  T1100 

0101 

T0001'  T0100'  T1001'  T1100 

1100 

T1000'  T0100'  T1010'  T1001 

0011 

T0001'  T0010'  T1010'  T1001 

0110 

T0010'  T0100'  T1100'  T1010 

1001 

T1000'  T0001'  T1100'  T1010 

1111 

1 

1 

T1000'  T01 00'  T0010'  T0001'  T1100'  T1010'  T1001 
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3  =  min[C^,C2,C3,C^,C1  2'C4  5^ 


and  branch  bi(Sg,S^,S2) 


dominates  branch  b2 (Sq,S3 ,S^) 
if  C1  3  t  p2  *  niin  [C^  ^C2  2^  P4  *  min[C4,C 

+  p^  •  min [C^ ,C3 ,C^  4,C3  2]  +  P3  *  min [C^ ,C3 ,C^  4»C 

—  2  +  P2  ’  min lci fC2 /Cf  4]  +  p^  •  min[C^,C2 

+  p3  *  min[C^,C3,C^  3,C^  4^  +  P4  *  m^n  3' 

or  Cl,3  +  P2  ‘  min  tc4  *'ci  r  2 1  +  P4  ’  min[c2*ci,2^  +  P 

+  p3  *  min  [Cx  ,C1  2J  <  C-^  +  p2  *  min  ,  3]  +  Pj^ 

+  p3  •  min[C4,C1^3]  +  p4  •  min[C3,C1^3l 

or  [Clf3  <  cost  <  C12]  <  Cl2 


2'C1,4'C1,2^ 

1,2] 

'C1,3'C1,4] 

C1,4J 

1  •min[C3,C1^2l 
•  min[C2,C1^3] 


which  is  satisfied  only  because  3  =  min  [C^] 

ieT(SQ) 

which  implies  that  this  is  the  only  condition  required  to 
guarantee  that  branch  b^(Sg,S^,S2)  dominates  branch  b2 (Sq,S3,S4) . 

By  the  same  procedure,  it  could  be  proved  that  branch 
b^(SQ,S^,S2)  dominates  also  branch  b3  (Sq,S,- ,Sg)  only  if 

Cl,3  =  min  t5i] 

ieT (SQ) 

Branch  b^(SQ,S^,S2)  dominates  br  ch  ;Sg,S3,Sg)  if 
C1  3  +  P2  ‘  min[C4,C1>4]  +  p4  *  minfC^C^]  +  Pj^  •  minfC^C^] 

+  p3  •  minlC^Cj.  41  <  +  p2  *  minlC^C^]  +  p3  •  mintC^C^] 

+  px  •  min[C4,C1  ,]  +  p4  •  min  iCj^Cj.  3J 
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or  if  IC^  3  <  cost  i  ci  4I  <  ci  4 

which  is  satisfied  only  because  ^  <  cost  of  any  other  test. 

Therefore,  the  branch  which  is  generated  by  test 
T1010  where  T^q^q  e  t  which  has  the  minimum  cost  at  node  N(Sq) 
dominates  any  other  branch  generated  by  tests  belong  also  to 
set  x. 

Eventhough  all  the  dominance  rules  are  proved  in 
case  of  having  five  or  six  LRUs,  they  could  be  considered  rea¬ 
sonably  as  general  cases.  However,  because  of  the  dramatic 
increase  in  the  number  of  possible  tests  (2^n“^  -  l)  in  case 
of  having  more  than  six  LRUs,  the  proofs  in  these  cases  are 


omitted  here. 


on  n 


c 


APPENDIX  d 


C  ********************************** 


C  MAIN  PRCoRAM 


C  ********************************* 


C  PU^PLii 

C  T  rf  1  6  PROGRAM  USES  A  BRANCH  AND  BOUND  ALGORITHM  TO  UNO  THE  OPTICAL 
C  SEOUtNCu  uK  TFSTS  FtGUIkEO  TO  LOCATE  A  MAuFUNCT  1  GNEO  UNIT  IN  A 
C  SYSTEM  UE  N  LFU'S 


INPUT 

CuNTRDL  CARDS 
DATA  CAkDS 


C  CUNTkOL  CAROS  (A) 
C  COLUMNS 


C  1-5  .....N . NUM3LR  OF  LRU‘S 

C  b-10  ...  .  .M  .  .  ..  .NUFQCk  OF  TESTS 


C  CONTROL  CARDS  <H)...(QNE  PEN  EACH  TCST) 

C  1-10  1 XX<  1 . I  * J ) . .TEST  1  DESCRIBED  BY  THE  N  BITS  (tACH  HIT  15 

C  REPRESENTED  3Y  J  ) 


C  DATA  CARDS  (AI...7F10.7 

C  COLUMNS 

C  1-10  . . bi'7C...C(l)  ..COST  OF  TEST  I 

C  DATA  CARDS  (H)...7Fl0.7 

C  COLUMNS 


C  l - l 0  . . . .P(J).., PROBABILITY  CF  FAILURE  OF  LRU  J 


C  DATA  CARDS  IO...7I10 
C  COLUMNS 

C  . . .... .61-70. . .MH(J) ...SECONDARY  ISOLATION  COST  Ct-  LRG  J 


C  PARAMETERS 
C  N=NUMBER  OF  L.R.U. 

C  M=NUMBER;OF  TESTS 

C  PI  J  )  =  PKU6AB1L  ity  of  failure  ur  L.R.U.  j 
C  Mm  J)  =  SECONDARY  ISOLATION  COST 
C  C< 1 J=COST  OF  TEST  I 
C  IS=  NODE  NUMBER 
C  lY(lS.J) =STAT E  UF  NUDE  IS 

C  1XX(1S,I,J)=  TEST  OF  MINIMUM  COST  AT  BRANCH  I  FROM  N^Dt  IS 
C  IXX1  l  IS  .1  •  J)  COMPLEMENT  OF  TEST  IXX(iS.l.J) 

C  CA(  IS.  I  )=C  OST  OF  TEST  OF  MINIMUM  COST  AT  BRANCH  I  F  REV  NUDE  IS 
C  LtIS)=LEvEL  OF  NODE  IS 
C  1 SS! I S) =PREVI OUS  NODE  TO  NODE  IS 

C  MIS):  NUMBER  LF  UNTESTED  L.R.U.  *S  AT  NODE  IS 
C  ALEI IS)=LO*tR  BOUND  AT  NODE  IS 
C  I D (  IS) . IF  IDI  I  S  )  =  1  ...NODE  IS  IS  DUMMY 

C  INDEX (  IS  *  I  )= 1  . BRANCH  I  FROM  NODE  IS  COULD  BE  USED 

C  INDEX!  IS, I  )  =  0  . BRANCH  1  FROM  NUDE  IS  COULD  NdT  t>c.  USED 

C  F I  IS.  I  )=VALUE  OF  FIGURE  OF  MERIT  TO  BE  USED  INRRANCHING  RULES  FOR 

C  BRANCH  1  AT  NODE  IS 

C  dUHDI  IS)  =  1  •  •  •  I F  THE  DUMMY  NODE  IS  IS  UPDATlD.ZERO  LlHlXtiJaE 

Dl MLNS 1UN  LYS1  (9) . CRRI 255) .LYS2C  9 ) ,L(  12)  ,  JUPD I H ) 

COMMON/ ADC/ 1  XXII  1  ,255.9) 

COMMON /BBC /C A I  11.255) 

COMMON/ OFF/  IXX II  1  1  ,255  ,9  ) 

COMMt.N/LMN/1  V<  1  1  ,9) 

COMMON/XY 2/  I  SSI  1  1  ) 

COMMUN/JVa/FI 11,255) 

COMMON/oH I /R I  1  1 ) ,  IDI  1 1 ) 

COMMON/ ENT /N  I  )  1  ) 

Common/ JMK/PI9) 

COMMun/BAB/ INDEX! 11,255) 

CoMMON/NAL/AL  61  1  1  ) 

CUMM0N/MAG/CI255) 

CCMMuN/ ANZ/PR 19 ) 

COMMON/ SEC/MH  19 ) 

CUMMON/ANA/CW I  255) 

CALL  DAT  A  I N, M ) 

C  AT  THE  F IRST  NODE 
1S=  1 

NUMBERS 1 

nude  *1 

ISMAXs 1 
L! 1 S) =0 
K I  1  S  )  =  N 
1DUS  )  =  0 

C  FIND  UOBER  BOUND  UB 
DU  162  J= 1 , N 
162  I  Y  I  1  S,  J  )  =  1 
SE  CC  =  0,0 
Do  29  J=1 ,N 

29  SECC  =  SECLFPI J)*MH I J) 

SUMB  =  0  •  0 
DC  30  1=1. M 

LAI IS.I  )  =  C!  1) 

CRRI  1  )=  C(  I  ) 

I  TT=0 

DO  31  J  =  1  . N 

31  I T  T  =  I T  T  ♦  I  X  X  (  1  ,  I*»  J  ) 

IF  I  ITT. EC. 1 ) GO  TO  32 
CRRI 1  )=0.0 
SUMO=SUMB*CRRI  I) 


32 


30  CONTINUE 

CMAXI=CRP (1 ) 

DO  35  I=2,M 

1FCCMAX I-CRRI  I  )  133.35. 35 

33  CMAXJ=CRR(J) 

35  CONTINUE 

S  UMCC=  S  UK  B  -  CM  AX  I 
UB  =  SUMC  C  + SECC 
C  FIND  L0*LR  BOUND 

CALL  REXP(IS.N.M) 

ALB (  IS)=R(  ISH-SECC 

C  STOPPING  TEST  USING  SECONDARY  ISOLATION 
H  SUM=  0, 

DO  37  J= 1 . N 
37  HSUM=HSUM+MH( J) 

IF((R(IS)*SECC).CT.HSUM)  GO  TO  34 

UB=HSUM 

WRI  TE  (6  »36  I 

36  FORMAT! /// »5X»  'USE  THE  SECONDARY  ISOLATION  F  OF  ALL  L.R.U.'I 
GO  TO  2752 

34  to RI TE (  6  » 85  9  ) I S»  ALB (IS) 

859  FORMA!  (  //,  I  OX  ,  •  I  S=*  •  I  3. 2X  .  'ALBUSJs'  ,F1  9.3) 

WR I TL ( 6  .222  >UB 

222  FORMAT  (  //  ,  1  OX  .  •  UB  =  ’  »F  10*3) 

1003  CALL  DOMINA(IS.N) 

CALL  FIND F ( IS.N) 

2003  CALL  BRANCH!  IS  »  ICP1  > 

WRI TE(6 .380 ) 

383  FORMAT ( /// *  i 0 X . 1 l S* . 20X .  •  NODE  •  .25X, • TEST*  ) 

WRITE!  6 .381 ) IS, ( < IY( IS. J) ,3= l ,N) . ( IXX< I S.1GPT , j). J= 1 ,N) ) 

381  F  ORMAT ( 1  OX, 1 1  ,21 X  ,9  I  1  .20 X  ,91 1  » 

C  FIND  T He.  Two  POSSIBLE  STATtS  FROM  THE  MOST  PROMISING  ERANCH 
MR  =  1S 
KS1  =  3 
KS2  -0 

DO  1  J=1.N 

LYS1(J)  =  IY(IS.J)*IXX(  IS. I  OPT, J ) 

LY  S2 1 J )  =  1 Y (  IS, J )*IXX l ( IS . I OPT .  J) 

IF(LYS1 ( J)-l) 9.4,9 
4  K  S I  =KS 1 ♦! 

9  IF(LYS2( J) .EQ.O)  GO  TO  I 

K52=KS2* 1 
1  Cunt  I nj E 

1  S=I  SMAXF  1 

IF (RSI • EQ • 1  )  GO  TO  6 
IF ( KS2 • Ed • 1 )  GC  TC  8 
C  USE  DUMMY  node 

WR I T  E( 6  » 3  83  )  IS 

383  FORMAT (//» 5 X NODE*  .2X.14.2X,* IS  DUMMY*) 

I D(  IS  >  =  1 
JUPD( IS  >=  0 
L(  1  S)  — L  (  MR  )  Ml 
R(  I  5  )=  3  • 

I SS( IS) =MR 

C  FIND  THE  PARAMETERS  OF  THE  TWU  NODES  BRANCHED  FRCM  THE  DUMMY  NODE 
IF  (KS2.GT  .KS1  )  GO  TO  7 
1 5=1 SF1 
DO  533  J=1.N 
500  I Y ( IS, J)=LYS2 ( J) 

L< I S) =LCMR)+2 
M  IS  )  =  XS2 
I D ( IS)=0 
1SS( l S)=l S- 1 
CALL  REXP(IS.N.M) 

CALL  LOWER!)! I S— l , N, I OPT ) 

15= ISF1 
DO  600  J=1  ,  N 
603  IY(  IS,3)=LYSI<  J) 


MIS)  =K  SI 
ISS( I  S  )  =  I  5—2 
L(  IS)*L(MP)*2 
I0(  I S )  =  0 
NODE=NODE*3 
1 SMAX=I SMAX*3 
CALL  kEXP(IS.N,M) 

CALL  LOKERBI  IS.N.  ICPT  ) 

GO  TO -  BOO 
IS*IS*I 
00  3  J- 1  ,N 
1 V I  IS. J)=LYS1(  J) 

L  (  IS  )=L  (MR  )*2 
K(  I  S)  =K  SI 
ID! IS)* 3 
ISS(  I  S  )  =  1  5-1 
CALL  kEXPIlS.N.M) 

CALL  L0*ERE3(  IS-  l.N.  IOPT ) 

I  S  =  I S+l 
DU  40  0  J= 1 #  N 
>0  I  Y (  IS. 0  )  =  LYS2 ( J) 

M  I  S)=KS2 
ISS(  IS)* I  5-2 
L (  IS) =L  t  MR ) *2 
I O ( 1 S)=0 
N0DE=N0DE+3 
1 SMAX=I SMAX+3 
CALL  REXP(IS.N.M) 

CALL  LOkERBI IS.N, IOPT ) 

GO  TO  800 

f  IND  THt  PARAMETERS  OF  THE  NEK  NODE 
DO  700  J=l . N 
>0  1 Y(  IS, J)=LYSI( J) 

M IS) =KS 1 
GO  TO  10 
DO  900  J* 1  ,N 
>0  I Y  l  IS, J )=LYS2  (  J) 

K(  I  S)=KS2 
)  ISS(  IS )=MR 

L  (  I S ) =L ( MR )  +  1 
ID( IS)=0 

CALL  REXP( IS. N.M) 

CALL  LOWERUd  S.N,  IOPT  ) 

NUD£=N0L>E* 1 
I  SMAX-=  |  SM  AX  ♦  1 

STuPPINU  TEST  USING  SECONDARY  ISOLATION 
>0  SUMPMH*  0  • 

SOM  MH  =  0 . 

DU  140  J=1,N 

IF(  IY< IS, J  ).EO.O  )  GO  TC  140 
SUMPMH*  SUMPMH* P(  J ) »MH( J ) 

SUM.4H*  S'  >MM  H*M  Ml  J  ) 

» 0  CONTI 

ir ( k<  | S  .  5UMPMH- SUMMH )  16  5*  lft» 1 06 
.6  bL  fc  =  ALB  (  IS  )-R(  IS  )  -SUMPKH+SUMMH 
IF (QLb-OB J 400  0 , 1  Oh, I  68 
40  00  WRITE <6,976)1  I  Y(  IS,J).J  =  1,N) 

9  ft  FORM  ATI ///,  I  OX,  •  USE  SEClnDARY  ISOLATION  TOK  UNTESTED 
2  A  T  NUDE •  , 2  X  , 2  0  I  1  ) 

UK  IT  EC  6.  fc>/0  )bU) 

TO  F  CRM AT (// , 1  OX  ,  • bLU  =  *  ,2X .Ft  .3 ) 

M  I*  (  2**1  X  (  1  S)-l)  )-l 
DO  lc>9  1*1 , Ml 
s9  M  I  S  ,  1  )  =  0  • 

GO  TO  167 

bS  IFIALbl  IS)-U8)333 ,168.168 
a  3  CALL  GENukA(  I  S.N.M) 

1F(MIS)-2)1000.  1  72,  1000 


1  T 


C  LAST  N  ODE  IN  THE  BRANCH 

172  SUMP3=0. 

DO  173  J=1,N 

IF(  IYI  IS,  J  1.EQ.0  1  GO  TO  173 
SUMP3=SUMP3«-P(  J) 

173  CONTINUE 

C  FIND  THE  ACTUAL  EXPECTED  COST  OF  THIS  BRANCH 
BN=( SUMP3 ) *CA ( IS,  1) 

BL  B=  ALB  {  IS  )-R<  IS  )+BN 
*R1 TE(6 .380) 

WRITE ( 6,361  II  S, (  (  IY(  1 S*  J )  *  J= 1 »Nl  *  ( I  XX (I  Sal  • J )  , J=1  ,N) ) 

WK IT  E ( 6 ,608  )ULB 

605  FORMAT! //, 1 0X,*6LB=* .F 1 0 . 3 ) 

IF ( BLB-Ub  1  167,168,166 

C  GO  TO  OTHER  BRANCH  OF  THE  DUMMY  NCDE  TO  ADJUST  Tf-E  LCaER  BOUND 

167  IF(L( ISJ-1) 174,175,174 

174  IF ( ID( ISS( IS)  )  .NE.l )  GO  TO  176 
1F( 1D( I S-l  l.EG.l I  GO  TO  1 7o 

GO  TO  177 
178  1S=ISS(1S> 

GO  Tu  167 

176  IS=IS-1 
GO  TO  167 

177  I S= I S-2 
JUP0(  IS  )  =  1 

ALB ( I  SI =BLB-R(  IS*1  ) 

1S= IS*I 

CALL  LOwERBi IS ,N,  I  OPT  ) 

GU  TO  800 

C  FAl HUM  THIS  NODE 

168  Ml=(2**  (MlS)-in-l 
DU  180  Is  I .MI 

160  r<  IS. I  1*0  . 

MR1 TE (6,330) 

330  FORM AT ( // ,  1  OX ,  *F ATHOM  LAST  NGDc  AND  CANCEL  LAST  T  t  ST  •» 

IF ( ID( I SS( IS ) 1-1 )20 1 . 181 . 20 1 

606  IS* I SS<  1  SI 
GO  TO  170 

181  M 1 * ( 2*  * ( K(  1  S) -1 ) ) -1 

DO  182  1=1. Ml 

182  F ( I S-l , I) =0. 

IS=  IS-2 

GO  TU  l 70 

C  A  FEASIBLE  SOLUTION 

175  UB=ULd 
MklTL(6.222)UB 
WR 1 TE ( 6 , 392) NODE 

JS2  FORMAT (//, 3X, •  NUMBER  OF  ACTIVt  NUDES  IS'.iei 

NUMJLR=NUMBORF 1 
Io= ISMAX 
C  bACKTRACR 

170  IF  (L(  ISI.FO.O)  GO  TO  299 
IF  (L  (  IS  ).EG  .1  I  GO  TO  200 
IF (  ID( 1 S) ,EO. 1  )  GO  TO  201 
IF(K(  1 S  )  ,£  O  .2  )  GO  TO  201 
Ml =(2** (K ( IS J -l I I-I 
DL  202  1=1 .Ml 

IF(F(  IS.  I  )  .NE.O)  GO  TO  300 
202  CONTINUE 

IF (  ID(  1 SS(  IS)  ),NL.l)  GO  Tu  201 
IF  ( JUPD ( ISS (  IS  )  I- 1 1666,20 1. 666 
201  I  S  =  I S-l 

GO  TO  1  70 

200  If  (  1D(  1  S)  .EC.l  )  GO  TO  204 
1F(  K(  IS)  .E0.2  )  GO  TO  204 
M1  =  (2**(K(  ISl-l)  ) -  1 
DO  205  1=1  .  Ml 

IF(F <  IS. I  I  .NE.  0. )  GO  TU  300 


20  5  CONTINUE 
204  IS=1 

299  00  20t>  1  =  1  .M 

1F(M  I5.D.NE.3.)  GO  TO  300 
206  CONTINUE 

GO  TO  2752 

300  ISMAX=IS 
GO  TO  2000 

2752  WRITE.16,391  )UB 

391  FORMAT! /// ,5X ,* THE  OPTIMUM  EXPECTED  COOT  IS  *,F10.5) 

WRITE (6, 3 92  INODE 

STOP 

END 

C  *********************************** 


SUBROUTINE  DATA 


*************  ********** 


********** 


PURPOSE 

DATA  IS  USED  TO  READ  bOTH  CONTROL  AND  DATA  CARDS 
SUBROUTINE  CAT  A  (  N  »  M  ) 

COMMON/ ABC / I  XX (  1  1  .255.9) 

COMMON/DEF/ IXX  1(  1  1,255.9) 

COMMQN/ANZ/PR (9) 

COMMON/* AG/CI 255) 

CUMMON/JHA/P< 9  ) 

COMMGn/SEC/MH(  9) 

READ  (5. 1 9  0  )N»  * 

190  F  CRMAT (215) 

»RITt(6. 191 )N.M 

191  FORMAT (  1M1  .//, IX.  ‘NOMOtR  OF  L.R.U.  I S • •  I  2 ,  1 0 A ,  • N JMb L R  L.F  TESTS 
IIS  *.14) 

WRITE! 6. 301  ) 

3  01  FQh MAT  1 ///  ,  6X  ,  •  I  *  *10X.'JXX(1,  I.J)*) 

Du  304  I=1,M 

RE  AD! 5 ,  302)(IXX(  l.I.J  ) »  3  = l.N) 

302  F  ORMAT (1011) 

WRITE! 6. 977 ) I , ( I  XX!  1  *  I • J)  , J  =  1 .N> 

9/7  F  ORMAT !5X. 13.10A .1011  ) 

304  CONTlNUt 

F  1ND  IXX 1 (  1 ,  I , J  ) 

DO  160  1=1  .  * 

DO  Itl  J  =  1  •  N 

IF | IXX!  1 . I ,3)  .EG.l  )  GO  TO  163 
1 XXI ( 1 , 1 . J) =1 
GO  TO  161 

1 o3  IXX1 (1 . I ,J)=0 
161  CONTINUE 
ItoO  CONTINUE 

WR1 Tt(6 ,1 95  > 

195  FORMAT! 1M1 ,//,20X, ‘COST  OF  TESTS*.//) 

DO  194  1=1. M 

RE  AD (5,6000)C!  1 ) 

6000  FORMAT!  7F 1 0.7) 

WR1TE(6, 193) I.C!  1 ) 

193  FORMAT!  1  OX.  *CCST  OF  TEST  •  .14.2X,'  IS*  »2X,F  10.7,/) 

154  CONTINJE 

WR1 TE(6 .196) 

196  F  ORMAT!///,  fcX  ,  •  J  *  .  1  3X  ,  •  F'  *  ,  1  5  X  ,  '  MH  •  ) 

READ  (5,192)  !P(J).J  =  1.N) 


on 


I  192  FOHMAT  (  7F1  0.7  ) 

READ  (5*6001)  (MH( J) « J  =  1 ,N) 

60  0  1  F 0  RM  AT ( 7 1  1 0  ) 

■  00  1 97  J»l  ,N 

WR 1 TE( 6.1 99) J,P( J) *MH( J) 

1**  FCRMAT(5X,  I2.5X.F 10.7, 5X,  I  10) 

197  CONTINUE 
RETURN 
END  - 

c  ******»***************»***w****«*4.* 


C  SUBROUT  INF  L  U  t»  E  R  £> 


C  *********************************** 


I 


PURPOSE 

LUAERB  IS  USED  TO  FIND  THE  L 0 MCF  BOUND  AT  nUUL  IS 
SUBROUTINE  LO*ERB< IS. N, 1UPT > 

COMMUN/LMN/I Y( 1 1 ,9) 

COMMON/KAL  /AL  B(  1  1  > 

CoMMON/GHl/Rt II). I D( 1 1 ) 

CUMMoN/XYZ/l SS( 1 1 ) 

COMMON/ t>PC/CA(  1  1*  255) 

COMMON/ Jr»K./P(9  ) 

COMMUN/AnZ/PR ( 9) 

DO  A99  J=  1  «  N 
♦99  PK(J)=P(J> 

IS1=ISS( IS) 

DO  40  J=1 .N 

IF(  IY(  1S1*  Jl.ca.il  GO  To  43 
PR  (  J  )  =  0  • 

43  CONTINUE 

C  SUMP 1  =  SUM  OF  THE  PRUBIL1TY  OF  FAlLURt  OF  THE  UNTLSltO  L.R.U.  AT  THE 
C  PRtVIOUS  NUOc.  To  NODE  IS 
SUMP 1=0. 

DO  44  3=1. N 

44  sum pi =  sump i ♦ pr ( j ) 

IF  (  1 D (  IS1  )-l>  48,47,46 

47  C1SS=0. 

GO  To  45 

46  C I S  =  =  SUMP 1  *C  A  (  ISI.IOPT) 

IF (  1 D ( 1 S ) • EO  *  1  )  GO  TC  46 

45  ACB(  IS)=AoS(lSl)-R(ISl)4R(IS)+ClSS 
GO  Tu  4  9 

46  AL3( I S) =ACB( ISl  ) -R( I  SI  ) 4R ( IS) 4CISS*R( 1S41 ) 

49  ARITC (6.659)1 S,ALB( IS) 

65  9  FUhMAT (// , 1  OX . • I S  =  *  .  13.2X  ,  • ALB (  1  j )  =  *. FI  0 .3) 

RL  TURN 
END 

C  ••♦••••*****»*****4**4**»*4****»*** 


I 

•4 


SUUROUT  INF 


BRANCH 


I 


PURPOSE 

FIND  THE  BRANCH  IOPT  »ITH  MAXIMUM  VALUE  OF  F(IS(I) 
SUBROUTINE  BRANCH ( IS* I OPT  J 
COMMQN/UVfc/F  (  I  1,255  ) 

CUMMON/ENT/M  1  1  ) 

M1  =  (2**(M  IS)-1I  )  -  1 
FMAX=F( IS , I > 


DU  96  1=2, Ml 
IF  (  FMAX  -F  <  IS,  I  )  )  95 
FMAX=F  (  IS,  I  ) 
CONTINUE 
DO  94  1  =  1  ,  MI 

IF  IF  <  IS,  I  )  •£.  Q,F  MAX) 
CUNT  INJE 
I  OPT— l 
F  (  IS.  I  )  =  0 • 

RETURN 

END 

******444* 


95. 96,96 


GU  TO 


SUBROUTINE  FINOF 


**♦*  *  4**44***  4  *  4  *  4  4  *  44*  *4  44  4  44  44  ♦  4 


PURPOSE 

FINDF  IS  USED  TO  FIND  THE  VALJES  Ur 
AT  NUDE  IS  »HICH  WILL  BE  USED  IN  THE 
SUBROUTINE  FINDF(IS.N) 

COMMON/ ABC/ IXX ( l  1  ,255,9  ) 
CUMMON/BBC/C At  1 1  ,255 ) 

COMMON/U V */F ( 1 1 , 255 ) 

COMMUN/LMN/IY  (11  ,9) 

COMMON/JMK./P(  9) 

COMMON/ BAB/ INDEX (11, 255) 
COMMON/ENT/M  1  1  ) 

WRi  Tfclo,  863) 

>3  FORMAT ( /// » l l X . •  I  •  ,23X , •  F(lS.l)*) 

M|s(2**(K< I S)-l) )-» 

PP 1 =  SUM  OF  THE  PROU.  OF  FAILURE  OF 

PP2=SUM  Oc  THE  PROB.  CF  FAILURE  OF 

TEsT  IN  BRANCH  I  PASSES 

PP=PRGoABILITV  THAT  THE  TEST  IN  BRAN 


F(IS.I)  U)  ALL  LkANCHl 
BRANCHING  PULES 


PP2  =  SUM  Oc  The  PROB.  CF  FAILURE  CF  ALL  OF 

TEsT  IN  BRANCH  I  PASSES 

PP=PRGoABILITV  THAT  THE  TEST  IN  BRANCH  I  WILL 
PPI=0. 

DU  90  J  =  I , N 

IF  (  I  Y<  I  S.  J  l.t'0.0  )  GO  TO  90 
PP  l  =  PPl FP( J  ) 

)  CUNT  INJE 

DO  91  1:1,  Ml 

IF(  INDEX (  IS,  I  )  *E  Q  *  0 )  GO  TO  93 
PP2=0. 

DO  92  J  = 1 , N 

IXYY=IY(  I5.J)*IXX(1S,1.J> 

IF ( IXYY.EO, 0)  GO  TO  92 
PP2=PP2  FP ( J ) 

>  CONTINUE 
PP=PP2/PPl 

F(  IS,  Il  =  -(PP*ALOG{PP)/.69  3F(  l.-PP)*ALOG( 1. 
GO  TO  864 
i  F<  IS. I  )  =  0. 

>4  *H1TC(6,862)I»F(  IS,  I) 


UNTE  STtO 
UNTESTED 


PASS 


-PP)/.693)/CA(  I  S  ,1  ) 


862  FORMAT! 10X, I4,20X,F10.4> 

91  CONTINUE 
RETURN 
END 

C  *********************************** 


SUBROUTINE  REXP 


********* 


************************  * 


TO  4  1 


PROBABILITY 


FAILURE  AMONG  THE  UNTESTED 


PURPOSE 

RfcXP  IS  USED  TO  FIND  THE  MINIMUM  EXPECTED  CUST  Tl  FIND  MALFUNCTIONED 
L.K.U*  FROM  NODE  IS 

SUBROUTINE  REXP(IS.N.M) 

COMMGN/LMN/1 Y< 11,9) 

COMMON/GHi/HI 1 1  )  ,  101 1 1  ) 

CUMMON/tNT/M  I  1) 

COMMON/ JHK./PI  9  ) 

C  CM  MON/MAu /C ( 255 ) 

COMMUN/ AN  A/CR ( 25S ) 

COMMON/ AN2/PR 19) 

DO  41  J  =  l  , N 

PR  (J)=P< J  ) 

IF  (  IYI  IS,  J  I.EO.i  I  GO  TO  41 
PKI J)=l . 

I  CUNT INUE 

P'41  N1  =  THE  minimum  probability  of  failure  AMONG  the  untested 

L  ,k  .U  .  AT  NOJL  I S 

PM 1 N2  =  T  HE  SECOND  MIN.  PROBABILITY  OF  FAILURE  AMONG  THE  UNTESTED 
L.R.U.  AT  NODE  IS 
PMIN1=PR( l ) 

DO  14  J-2.N 

IF  (  Pkt  J)-PMINI)  13,13.14 
)  PMIN1=PRIJ) 

1  CONTINUE 

DO  I  6  J  =  I  ,  N 

IF(PR(J).EO.PMINl  )  GO  TO  1? 
i  CONTINUE 

r  pr(j>= l. 

P  M  l  N2  =  P  R  (  I  ) 

DO  23  J  =  2 . N 

IFIPRI JI-PM1N2  )  19,19,20 

4  PM1N2=PR(J) 

)  Continue 

DO  599  1=1  •  M 
1 9  CR  (  I  )  =C  <  I  ) 

JJ  =  M-M  IS  )*1 
DO  27  1=1, JJ 
CM AX=  CR (  1  ) 

DO  24  J  =2  *  M 

IF C CMAX-CRC J) >  23,24,24 
>  CMAX  =  CMJ) 

1  CONTINUE 

DU  25  K Z-  I  , M 

IF (CRIK2) .EG.CMAX )  GO  TO  26 
a  CONTINUE 
"j  CR  (K2  )=  0  , 

7  CONTINUE 

C  SU  4=  SUM  OF  THE  (K(IS)-l)  MINIMUM  COSTS  OF  TESTS 
CSUM  =0 . 

DC  28  I =1 . M 


23,24,24 


MINIMUM  costs 


TESTS 


28  CSUM=CSUM*CH(  I  ) 

R(  IS)  =  (PMI  N 1  ♦  PMI  N2)*CSUM 
MR ITE( 6  * 839  )  J  S  ,  R  ( IS) 

839  FORMAT  C //.  1  OX,  •  I  S=*  . 13 ,3X , • R <1 S )  =  • . FI  0 .3  ) 

15  RETURN 
END 

C  *********************************** 


C  SUBRUUT  I  NC  D  O  M  I  N  A 


C  *********************************** 


C  PURPOSE 

C  DOM IN A  IS  UStU  TO  FIND  1 NDE X (  I S  .  I  )  AT  NODE  IS  FLR  ALL  POSSIBLE 
C  BRANCHES 

SUBROUTINE  DOMINA(IS.N) 

C  IF  INDEX!  IS, I  1  =  0  ...BRANCH  I  IS  DOMINATED  BY  ANOTHER  ERANCH  WHICH 
C  ITS  INDEX  EQUALS  I 

0 1 ME NS I  ON  PT l ( 255)  , PT 2 1  2 55 ) , PT ( 255 )  ,K S3 1  ( 2 55  )  ,KSS2< 255 ) ,CAA( 255 ) 
DIMENSION  IV  SSI  l  R),  I Y  SS2(  9  )  ,  P  A  A<  9  ) 

COMMCN/AbC/I XX( 1 1  ,255.9) 

COMMON/ BBC /C A (  I  1 .255) 

CCMMON/DEF/ IX X  I (  1  1  ,255,9) 

CGMMON/LMN/I Y( 1 1 ,9) 

COMMON/ EnT/K (11) 

COMMON/ JHX/P(9 ) 

COMMON/BAB/ I NDE  X (  11,255) 

M1=12**(K ( IS)-1) )-l 
DO  50  I *1 .Ml 
PT1 ( I)=0. 

PT2  (  1  )=0. 

K SSI  (  1  )=0 
KSS2 ( I ) *0 
DO  51  J  =  1  ,  N 

1VSSUJ)=1Y(1S.J)*IXX(IS.I,J) 

I Y  SS2 ( J  )  -  I  Y (  IS.J)»1XXI(1S*1.J) 

IF  (  l  YSSU  J  ).FQ.  0)  00  TO  52 
KSS 1 i  I  )=KSSI(  1  )♦  1 
PT 1 ( 1 )=PT1 ( 1 ) ♦  P  (  J  ) 

52  IF  (  1YSS2(  J  ).E  J.O)  oO  To  51 
KSS21 1  )=N5S2<  I  )  ♦ l 

PT2( 1)=PT2( 1)*P( J) 

51  CONT  lNUE 

IF(KSS) (1 )-KSS2( I ))  54.53,55 

53  IF(PT1(  I  1-PT2U  )  )  5*. 5«*, 55 

54  PT ( I )  =PT 1(1) 

GO  TO  50 

55  PT (  I  >=PT2<  I  ) 

50  CONTINUE 

DO  60  1 =1 ,M1 

60  C AA (  I  )  =  CA (  IS,  I  > 

C  FIND  MINIMUM  COST 

CAAMIN=CAA<  I  ) 

DO  62  1=2, Ml 

1 F ( CAA (  I  )  -C  AA  Mi N ) 61  ,61,62 

61  CAAM1N=CAA(  1  ) 

62  CONTINUE 

OO  63  1 = 1 • M I 

IF <CAA(  I  )  . EU.CAAM  IN  )  GO  TO  64 

63  CONTINUE 

64  CA l 1  =  CAA(  I  ) 


CAA( I )  = 100000. 

1 A  As  1 

DO  65  I  =2, Ml 
IF  ! CAA I  1 1-  C22  ) 66, 66,05 
66  C22—CAA { 1 > 

65  CONTINUE 

DU  70  J=1,N 
PA A !J)=P!J) 

IF<  IV!  IS,  J  l.EQ.l  )  GO  TU  70 
PA A ( J  )  = 1 . 

70  CONTINUE 

PA AM  JN=PAA(  1  ) 

OU  73  J=2 , N 

IF (PAA( 0) -PAAMIN  >  72,72,73 
72  PA AM I N*  PA A !  3  ) 

7 J  CONTINUE 

C  CHECK.  IF  THE  BRANCH  WITH  MINIMUM  COST  OF  TEST  DOMINATES  ALL  OT  HEW 
C  BRANCHES  OK  NOT 

IF  CCA  11 -PA AM I  NFC  22 >74, 74 ,75 

74  DO  76  1=1, Ml 

7t>  INDEX!  IS,  I  )  =0 

GO  TO  6  0 

75  1F1KSS1 (  IAA1-KSS2! IAA)  )  77,78,  77 

77  DO  79  1 =1 ,M1 

IFCCAU-PTI  I A  A  *C  A  C  IS,  I  J)  161,61  ,82 

81  INDEX! IS, 1 >=D 
GO  TO  79 

82  INDEX! IS, I  1=1 

79  CONTINUE 

Go  TO  8  0 

76  DO  83  1=1 .Ml 

IF! KSSl  !  I  > -KSS2! 1  )>  64.85, 84 
64  IF  I CA1 1 -PT  C  I AA*CA ( IS.  I) ) ) 86. Bo,  67 

eo  INDEX (IS, I )=0 

GC  TO  83 

87  INDEX! IS. I  1=1 
GO  TO  83 

85  IF!X! 1S)>4 >88,89,86 

88  INDEX! IS. 1  1=1 
GU  Tu  83 

89  INDEX! IS. I  1=0 

ej  CONTINUE 

eo  INDEX!  IS.  IAA  )=  1 

WETUKN 
END 

C  ********************************  * 


C  SUBROUTINE  GENBKA 


C  4***********4F****F****F***»***** 


C  PURPOSE 

C  GENBHA  IS  US FD  TO  GENERATE  ALL  POSSIBLE  BRANCHES  FROM  NODE  IS  , 
C  1  HEN  FIND  THE  TEST  WITH  MIN  COST  IN  EACH  BRANCH 
SUBROUTINE  GENBKA!  1S.N.M1 

DIMENSION  1XYXI255.9)  ,  KS  I  !  255  )  .  MMM!  255  > 

COMMON/ABC /I  XX!  I  1  .255,9) 

COMMON/ BbC/CA!  1 1 , 255  ) 

COMMUN/OEF / I  XX 1 !  1 1  .255  ,9) 

COMMON/LMN/ I Y ( 11,9) 

C  CM MON/ MAo/C ! 2 55 ) 


t 


CUMMON/ENT /K (  1  I  ) 

DU  100  1  =  1. M 

KS 1(1 )=0 
DO  1 01  J=l  • N 

IX  YX(  1  •  J  )  =  l  Y<  I  St  J  )*  1  XX{  1  t  I  .  J  ) 
IF11XYXU,  Jl.tQ.O)  GC  TU  101 
KSK  l)=KSi  (I)*l 

101  CONTINUE 
100  CUNTINUE 

I JK=(K( I S)/2>  *2 

IF(  IJK.EG.KI IS  )  )  GO  TO  110 

DL  102  1—1  *  M 

1F(KS1( I).LT.(K( 1S)+1)/2)G0  TO  104 
DU  103  J= 1 * N 

1 X YX( 1 ♦  J  )  =  I  Y(  IS* J)*1XX1  C  1  *1  * J) 

103  CONTINUE 

104  DO  105  J  —  l  *N 

IF ( 1XYX< i * J) .EO. 1 )  GO  TO  108 

105  CUNTINUE 

12 1  MMMU)  =  0 
GU  TU  1  02 

106  DO  10O  J=1  •  N 

IF  (  1 X YX (  l.J)-lY(IS.J))  109*1 Oo  *  1 09 
106  CONTINUE 
uO  TO  121 

109  MM*4(  I  )=  1 

102  CONTINUE 
GC  TO  120 

110  DU  130  1=1. M 

IF  (K51  (  1  )-K<  I  S)/2)  124.129*122 

122  DO  123  J= 1 *N 

123  IXYX(1*J)-1Y(  IS. J  >*1 XX  1(  1  *  I  * J ) 

124  DO  125  J=1  .N 

IF(  IXY  X (  I »  J  )  *EO  •  1  )GU  TO  128 

125  CUNTINUE 
131  MMMI 1 )- 0 

GU  TU  130 

128  DC  126  J=l  *  N 

IF  (  I  X  YX  (  I*  J  )-  1Y<  15. J))l  29  *126*129 

126  CONTINUE 
GU  TU  131 

129  MMMU)sl 

130  CONTINUE 
GC  TU  1200 

120  »H1TE(6.650) 

85 J  F  URMA  T ( ✓// .  1  1  X  *•  I*  *2  0X.«  IXXC  IS  •  1  •  J)  a 1 
1=  l 

DU  111  11  =  1  •  M 

IF  (FMM(  i  I  )  *C  i.O)  GU  TO  111 
K  Ms  1  I 
LLL= I  1*1 
DU  112  l=LLL*M 
IF ( MMM(L ) • EC.O )GC  TO  1 
DO  115  J=1 *N 

IF  C  I X Y X ( II  *  J  )  - IXYX(L  *  J  )  ) 1 12* 115*112 

115  CU NT  I NJ  E 
MMM( L) =0 

IF CCINM  |-C(L  )  >118*118*  116 

116  KM=l 

118  CA( IS. 1 )=C ( KM> 

DC  11V  J= 1  *  N 

I  XX ( I S. I  . J) =1 XX(  l  .KM. J) 

IF  (  IXXC  I  S.  I.J)  •EU.DGU  TU  5 
1XX1  ( IS* I  *  J  1-1 
GU  TU  119 

5  1 XX 1 (  IS*  I  *  J  )=  0 

119  CUNTINUE 

II*:  CONTINUE 


1  =  1*1 

in  continue 

GO  TO  1603 
1230  HR lTt( 6, 853) 

1  =  1 

DU  1111  11=1 » M 

1F(MMM(  II )  .EU.3)  GC  TO  1111 

KM=  I  I 

LLl.sU  +  1 

DO  1112  L  =  LLL t  M 

IF ( MMM(L) .EC.O  )GO  TO  1112 

DO  1115  J=1.N 

IF  <  IXYX<  II*J)-1XYX(L«J))  1  1 13»1115>l  113 

1115  CONTINUE 

1117  MMM(L  )=  3 

IF (C(XM  i-C(L)  )  1 1  ie , 1  1  18, 1  116 

1116  KM-L 

1118  C A (  IS.  I  )  =  C(KM  ) 

GO  TO  1700 

1113  DO  1114  ju=1 t  N 

IF ( IXYX (11*  JU )*IXYX(L*JU)“1Y(  1 5 1 JO ) )1  112*1 114*111 

1114  CONTINUE 
GO  TO  1117 

1700  DO  1119  J=1 .N 

1XX(  15* I  • J)=I  XX(  1  .KM, J) 

IF ( IXXC  IS.  1 . J  )  .EO.l )GO  TC  1035 
1XX1  (IS.l  , J)=l 
GU  TO  1119 
1035  IXX1 ( IS . 1 .J)=0 

1119  CONTINUE 
1112  CONTlNJE 

1=1*1 

1111  CONTINUE 
1633  MI  =  (2**(K(1S>-1>  )-l 
DO  652  1=1 i Ml 

I TE( 6  »  85 1  )  1  .  <  I  XX  (  I S. I • J  )  • Js|  ,N) 

651  F OKMAT< 1  OX  ,  14  ,20X  .9  1  1  > 

852  CONTINUE 
RETURN 
END 


